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An expression for the work function of a monovalent 
metal, previously obtained by E. Wigner and the author, 
included a term which represents the energy required to 
move an electron through an electrostatic double layer at 
the surface of the metal. In order to determine the moment 
of this double laver, it is necessary to make an explicit 
calculation of the electronic charge density at the surface of 
the metal, which is attempted in the present paper. In the 
model taken for this calculation, the density of positive 
electricity is assumed constant on the negative side of the 
YZ plane, and is zero in the positive half-space. The 
electrons are so distributed as to neutralize the positive 
charge density in the interior, and the metal as a whole is 


neutral. An approximate self-consistent solution of the 


Fock equations gives the charge density at the surface. An 


1. INTRODUCTION 
HE theory of the work functions of mono- 
valent metals has been discussed in a recent 
paper! by E. Wigner and the author. This work 
was incomplete in two respects. In the first place, 
the free-electron values were used for the Fermi, 
Coulomb, exchange and correlation energies of 
the This defect 
making explicit calculations of these energies for 


metal. could be removed by 


the different metals. Such calculations have been 
carried out, at least in part, for Na,?:* Li,‘ ° 


*Society of Fellows. The major part of this work was 
completed while the author was a Fellow in Mathematics 
at Princeton University. 
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ittempt 1s made to go bevond these equations, so as to 
include the effect of correlation (or polarization) forces on 
the double layer. The density of positive electricity assumed 
is the same as that which would exist in Na if the ions were 

continuous charge The final 
values obtained are 2.35 and 2.0 ev for the work function, 
and 0.4 and 1.0 ev 


where the first values in each are obtained when the correla- 


treated as a distributior 


for the moment of the double laver, 


tion forces are included, and the second when these forces 


are omitted in the calculation. A comparison of these values 


with the experimental values of the work function of Na is 
given. It is concluded that the surface barrier is due pri- 


marily to exchange and polarization forces, and that 


ordinary electrostatic forces play a minor role 


and Cu.® ? Secondly, the electrostatic double 
layer at the surface of the metal was left unde- 
termined. 

The main object of the present work is to show 
how the double layer is determined and to make 
a rough calculation of the moment of the double 
layer at the surface of an ideal metal. The charge 
distribution at the surface is calculated by the 
method of the self-consistent field based on the 
Fock equations. This calculation may be of some 
interest, not only in the present connection, but 
also because it yields some information about the 
potential barrier at the surface of the metal. 

In the previous paper, the work function was 


obtained by first calculating the energy E(n;, n.) 


®°H. M. Krutter, Phys. Rev. 
7K. Fuchs, Proc. Roy. Sox 


48, 664 
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of a metal containing m; ions and n, electrons. 
This calculation required but a slight modifica- 
tion of the method used by Wigner and Seitz® 
for calculating the energy of an uncharged metai. 
The work function, ¢, was then easily obtained, 
as 


g= —(0E(n;, n.)/ One) nan. (1 


It was shown that the work function consists of 
two parts, the “binding energy”’ of the electron, 
and the energy required to move the electron 
through an electrostatic double layer. at the 
surface of the metal. The first part may be cal- 
to the surface, and 
(I). 


determine the moment of the double layer, it is 


culated without reference 


this has been carried out in In order to 
necessary to make an explicit calculation of the 
charge density at the surface. 

As our model for the present calculation, we 
suppose that the metal occupies the half-space 
Y-Z plane. The 


density of positive electricity is constant for 


on the negative side of the 


x <0 and is zero for x >0. This same model has 
been treated by Frenkel’ from the standpoint of 
the Thomas-Fermi theory. We here determine 
the charge distribution at the surface from an 
approximate solution of the Fock equations. 
Our calculation also gives the work function of 
the model from a somewhat different viewpoint 
than that used in (I). The density of electrons 
assumed is approximately the same as the density 
of free electrons in Na, and a comparison is made 
of the work functions of our model and that of 
Na. An attempt is made to include the corre- 
iation (or polarization) forces which are neglected 
in the Slater-Fock approximation of one-electron 
wave functions. 

2. DEFINITIONS 


The work function may be defined (cf. I) as 
the difference in energy between a lattice with 
an equal number of ions and electrons, and the 
lattice with the same number of ions, but with 
one electron is assumed in both 
cases that the lowest electronic states are com- 
pletely filled, so that the electron is removed 
from the highest energy state of the neutral 


removed. It 


* Cf. reference 2. A general review has been given by J. 
C. Slater, Rev. Mod. Phys. 6, 209 (1934). 

* J. Frenkel, Zeits. f. Physik 51, 232 (1928). This paper 
also gives a discussion of the surface double laver. 
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metal. Since the work function may be different 
for different crystallographic planes, it is neces- 
sary to specify the position of the electron after 
its removal from the lattice. It is supposed that 
the electron is removed to a point in the neigh- 
borhood of a surface plane of the crystal, the 
distance from this plane being small compared 
with the dimensions of the plane, but large in 
comparison with atomic dimensions. The poten- 
tial differences between the outer neighborhoods 
of differently oriented crystal surfaces will be 


equal to the differences of the work functions of 


these surfaces. 

The moment of the electric double layer on 
the surface may be defined as the difference in 
potential between a point outside the surface 
and a point in the interior. This definition will 
serve for our free-electron model. As the potential 
in the interior of an actual metal is not constant 
in the interior, but is periodic, it is necessary to 
specify further the position of the interior point 
relative to the lattice. The position of this point 
is somewhat arbitrary, but once chosen, it fixes 
the otherwise arbitrary constant in the energy 
parameters of the wave equations for the dif- 
ferent electrons which are valid in the interior 
of the metal. It will be sufficiently accurate to 
specify the position of the interior point with 
respect to the approximate model of Wigner and 
Seitz in which the actual metal is replaced by a 
set of close-packed spheres (the s spheres).? In 
each sphere, the potential will be a function only 
of the distance from the center of the sphere. As 
our interior point we shall take, as in (I), a point 
on the surface of the s sphere. This definition is 
most convenient, as it gives a zero moment if the 
electron-distribution in the s spheres of the sur- 
face ions is the same as that in the s spheres of 
the interior ions. The double layers are due to 
the fact that the actual distribution may extend 
outside the limits of the s spheres of the surface 
ions, and may have greater or smaller densities in 
the interior of these spheres. 

The exact definition of the moment of 
double layer is of importance in our connection 
when we compare the work function and double 
layer of our free-electron model with the corre- 
sponding quantities of the sodium crystal. We 
shall find that the work functions may be 
directly compared, but that, as a consequence of 


the 
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the above definitions 
aver of the 110-face of the 


the moment of the double 
sodium crystal is 
about 0.25 ev less than that of the free-electron 
model. 

DouBLE LAYERS 


3. MAGNITUDE OF THI 


Before proceeding with the calculation, it 
might be well to say a few words about the ex- 
pected order of magnitude of the moment of the 
double layers. The most direct evidence was 
It was found that the theoretical 
values of the work the alkalies 


checked the experimental values very closely, 


given in (I 
functions of 


even though the double layer was omitted en- 
tirely. This suggests that the double layers of 
these elements are probably small (less than half 
a volt). 

Since the potential in the interior of a metal is 
constant (or, more accurately, periodic), the 
the work of different 
surface planes of the same crystal are equal to 
the 


differences in functions 


the differences in the moments of corre- 
sponding double layers. These differences, experi- 
mentally, are of the order of 3 to 1 volt.’ One 
might expect that the moments themselves are 
of this same order of magnitude. 

On the other hand, experiments show that the 
value of the work function depends to a large 
extent the | Ad- 


sorbed gas or deposits of foreign material may 


on condition of the surface.' 
alter the work function considerably. Ives and 
Olpin® and J. J. Brady’ have shown that a 
monatomic layer of alkali atoms deposited on a 
metal may reduce its work function by as much 
as 50-75 percent. These data indicate that the 
work function of a surface is really determined 
by the top few layers of atoms, and not by the 
metal as a whole. The large changes in the work 
function due to conditions, 
must, of course, be due to corresponding changes 


different surface 
in the surface double layer. A knowledge of the 
charge distribution at the surface is thus really 
necessary for the the work 
function. One may expect, however, that the dis- 


determination of 


Farnsworth and B. A. Rose, Proc. Nat. Acad. 


N. E. 
Sci. 19, 777 (1933); B. A. Rose, Phys. Rev. 44, 585 (1933); 
N. Underwood, Phys. Rev. 47, 502 (1935). 


'See, for example, A. L. Hughes and L. A. DuBridge, 
Photoelectric Phenomena (McGraw-Hill, 1932), p. 72. 
*H. E. Ives and A. R. Olpin, Phys. Rev. 34, 117 

J. J. Brady, Phys. Rev. 41, 613 (1932). 


1929). 


WORK FUNCTION 


tribution for pure metals will be such that the 
double layer is small. 
4. THe FREE-ELECTRON Mone! 

The density of positive electricity has the con- 
stant value ep, in the half-space on the negative 
side of the YZ plane and is zero elsewhere. If the 
electrons is the electrostati 


density of p(x 


potential, I’, will be given by the solution of 


—drep(x); p(x 


If the charge distribution is known, this equation 


may be solved for I’. The boundary conditions 


are, if the metal is uncharged: 


V=0, dV/dx=0 as 


dV /dx=0 as x—- 


Condition (3b) requires that!* 


lm p(x dx =. 4) 


— 
. 


The constant value of the potential in the in- 
terior,!> V°, is, by definition, the moment of the 


surface double layer. 


5. THE Fock EQUATIONS 


Our calculation of the charge density at the 
surface is based on the Slater-Fock'® approxima- 
tion of one-electron wave functions. An attempt 
is made, however, to correct for the correlation 
or polarization) forces which are neglected in 


this approximation (cf. section 8). If there are 


? 


electrons which occupy n 2 doubly degenerate 


states, the Fock equations have the form: 


(— (A? /2u)A;—eV(r;) —A,)¥,(r;) = EW,(r;), 5) 


where |’(r;) is the Coulomb potential of the com- 
bined positive and negative charges. If 


e(r, r;) = Dv* (ry i(r)) 6) 


‘ If the metal has a surface charge, the integral of Eq. (4) 
is no longer zero. In this case, condition (3a) is replaced by 
dV/dx=4n0 as xx, where oa is the surface density of 
charge. 

‘8 The superscript zero will be used throughout to indi- 
cate values which apply in the interior. 

6 |. C. Slater, Phys. Rev. 34, 1293 (1929); 35, 210 (1930). 
V. Fock, Zeits. f. Physik 61, 126 (1930). An excellent re- 
view has been given by L. Brillouin, Les Champs ‘‘self 
consistents”’ de Hartree et de Fock, Actualités scientifiques 
Herman, Paris, No. 159, 1934). 
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is the Dirac density matrix, the exchange oper- 


ator A; is defined by 


The sum in (6) is over all orbital wave functions, 
the spin factors being omitted. The factor 3 
appears in (7) because the exchange operator is 
effective only for electrons of parallel spin, while 
the wave functions of all electrons are summed in 
p(r, ¥;). The density of electrons is p_(r)=p(r, r). 

In the following, we shall use instead of the 
exchange operator, an exchange potential defined 


er 
A ,(r;)= | 
2 


This exchange potential thus depends on the 


by :” 
p(t, r;)¥;(r)dr 


r—r;|~;(f;) 


wave function, and will be different for different 
wave functions. The advantages of using the 
exchange potential rather than the exchange 
operator are two. First, the potential is easier to 
visualize physically. As Slater and Krutter!’ 
point out, the field acting on the electron is 
actually to be identified with the sum of the 
Coulomb and exchange potentials, rather than 
the Coulomb potential alone. There is a hole cut 
out of the charge distribution about each electron 
due to the fact that the probability of two elec- 
trons of parallel spin being close together is 
small. The total charge cut out by the hole is 
exactly e. The exchange potential is essentially 
the potential resulting from the hole.’ The 
second reason is more practical. We want to find 
a self-consistent solution of these equations, and 
this must be obtained by numerical calculations. 
The exchange integrals are very difficult to 
evaluate directly. The procedure adopted was to 
approximate the actual set of wave functions by 


‘7 Cf. reference 8, and J. C. Slater and H. M. Krutter, 
Phys. Rev. 47, 559 (1935). 

18 Tt should be emphasized that this interpretation can- 
not be made in all cases. The exchange potential defined 
in this manner may be complex or may have infinite 
values. In our case, the potential is real and finite. The 
probability of an electron in the jth state being in the 
volume element dr, and a second electron being in the 
volume dr, is the real part of: 


; ,* 


Di(Ws* (i We(ti dy ;* (nv, (8) —vi* (ty (ey; * (ev (01) drdr,. 


The second term represents the ‘‘hole.”’ 
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Fic. 1. Exchange potential in interior as a function of 1 
velocity of the electron. 


a set of analytic functions for which some of the 
required integrations could be performed. It was 
found that the exchange potentials were much 
less sensitive to the actual form of the wave 
functions than the exchange integrals. 

THE INTERIOR 


6. WAVE FUNCTIONS IN 


We shall replace the actual infinity of wave 
functions for the free-electron model by a finite 
number by introducing periodic boundary con- 
ditions (period L) in the Y and Z directions, and 
a fixed boundary at x= —L. The wave functions 
in the interior!’ will be the plane waves 


y, 3)=2°L 


Wi. (x, Psin (Rixty)e Atk), (9 
The allowed values of ko and k; are 2rve/L and 
2rv3/L where vo and pv; are integral. The allowed 
values of &,; are 2rv,'L where the values of v 
lie close to the half-integers (4, 1, --+), We 
suppose that there are » electrons in the volume 
L’ and that all states are doubly occupied from 
to where and 
Vm = (3n/8r)'. The radius of the sphere (s sphere) 


v=(0 V=Vm, v= (y+ vo? + v3”)? 


whose volume is equal to the average volume 
occupied per electron is: 
(10 


r,=(3/4an)'L =1.92(2zrv,,/L) =1.92R». 














THEORY OF THI! 
We shall find it convenient to use r, as a param- 
eter which gives, indirectly, the density of the 
electrons. 

The 
constant. The exchange 
wave function of the form (9) has the effect of 


the interior is 


operator acting on a 


exchange potential in 


multiplying the wave function by a constant. 
This constant, evaluated by 
Dirac,’ depends on the ratio a=k/k,, of the 


which has been 
momentum of the given electron to the maximum 
momentum of the Fermi distribution. The value 


of the exchange potential in the interior is: 
Ao =eknf(a)/2r=0.306e f(a) / rs, 


where 


f(a) =24+ (1/a)(1— a’) log ((1+a)/|1—aj). (11) 


A plot of A,° as a function of a is given in Fig. 1. 


° is: 


The equation for ¥ 
(— (h?/2u)A—eV°—A,)pP =Ew,’, 


so that E,.= (h?/2u)k?—eV°—A;”. 


It may be noted here that, to the approxima- 
tion of the Fock equations, the energy required 
to remove the &th electron from the metal is 
—fF,. The work function, ¢, is the minimum 


value of —£,, or 


g= —(h?/2u)kn2+eV?+Ax,”. (14) 


This agrees, except for the correlation energy 
which we have so far neglected, with the results 


of (I). 


7. METHOD OF CALCULATION 


The exchange potential decreases in magnitude 
as one approaches the surface, and goes to zero 
as x. It is independent of y and sz. The Fock 
equation (cf. Eq. (5)) thus has the form: 


(— (hk? /2p)A—eV (x) —Ax(x) )vi(r) =Epi(r) (15) 


with V(x) given by (2) and (3) and E, by (13). 
Since the potential is independent of y and z, 
the wave functions must be of the form: 


Vi(x, ¥, 3) = 9(R, xe truths? (16) 
The equation for ¢(k, x) is 
(— (h?/ 2m) (d?/ dx?) —e V(x) —Ax(x)) o(k, x) 

= (E,.— (h?/2u) (ko? +k3")) o(k, x) (17) 





’?P. A. M. Dirac, Proc. Camb. Phil. Soc. 26, 376 (1930). 
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and the asymptotic form of ¢g(k, x) in the in- 


terior is: 


o(k, x) =2'L-4 sin (Rk, x+y). 18 


The phase, y, may depend on ke and &; as well 


as k,. 
In order to an accurate solution of 
takes a 


wave functions, and then computes the Coulomb 


obtain 
Fock’s equations, one set of starting 
and exchange potentials from these. The next 
approximation to the correct wave functions is 
obtained by solving Fock’s equations using these 
potentials. This procedure is repeated until a 
self-consistent set of wave functions is obtained, 
i.e., a set such that the Coulomb and exchange 
potentials computed from these wave functions 
are the same as the potentials used in deriving 
the wave functions. Because of the numerical 
difficulties involved in evaluating the exchange 
integrals, it was not found feasible to carry out 
this process in full. 

Instead, the exchange potentials were chosen 
at the beginning and held fixed throughout. A 
self-consistent solution was then obtained for the 
electrostatic part of the problem. In the final 
stage, the Coulomb potential calculated from 
the final wave functions was equal to the poten- 
tial used in determining these wave functions. 
The exchange potentials were determined from 
an approximate set of wave functions which will 
be described in the next section. As the charge 
distribution corresponding to these approximate 
wave functions is fairly close to the calculated 
charge distribution, one may expect that the 
exchange potentials as determined from the two 
sets of wave functions are also very nearly the 
same. The Fermi hole is very large near the 
surface, and so should not be sensitive to the 
exact form of the wave functions at the surface. 
POTENTIALS 


8. CALCULATION OF EXCHANGE 


As one approaches the surface of the metal, the 
Fermi hole will increase in volume in such a way 
that a total charge e is always cut out by the hole. 
The shape of the hole will depend not only on the 
speed of the electron, but also on its direction. 
The dependence on direction is rather small, and 
is unimportant for our purpose. The exchange 
potential of an electron whose motion is per- 
pendicular to the surface will drop a little more 
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slowly than the potential of an electron of the 


same energy whose motion is parallel to the 
surface, but the drop will start sooner, so that 
the potentials of the two electrons will be about 
equal half-way down the slope. We shall, in our 
work, neglect this dependence on direction. 

It is very difficult to evaluate the exchange 
potentials for wave functions other than the 
plane waves discussed above. It is possible, 
however, to get some idea of the behavior of the 
exchange potentials by calculating A,(x) for the 
wave functions of a somewhat simplified model. 
We shall suppose that the electrons move in a 
constant negative potential, W, for x <a and that 
the potential is zero for x >a. The position of the 
barrier, a, and the height of the barrier, Il’, are 
adjusted in such a way that the wave functions 
and charge distribution of this simplified model 
correspond most closely with the wave functions 
of the free-electron model. The integrals are so 
complicated, even for this simple model, that it 
to the 


potential wall, x=a. 


is difficult evaluate them except at 
Instead of using x as a coordinate, we shall use 
&=x—da, so that the barrier is at £=0. The wave 


£ 
functions are then of the form (16) with 


2'L 1 p)(—(p —k’)* sin ké 
o(k, t)=4 +kcoskt) §<0, (19) 
L2IL-¥(k /p)e-(ez- kt z>0, 
where r= — Qu Wh’. (20) 


The integrals we must evaluate are those of 


Eq. (/). 


e- y *(r)Wi(r)ve(ri)dr, 
A.(i)y(r)=—> —- . (21) 
2 r—Tr; 


The summation over / is over all occupied states, 
and may be replaced by an equivalent integra- 
tion over asphere S of radius ,, = (27, L)(3n 82r)'. 
(n is the number of electrons in the volume L*.) 
It is shown in the appendix that if we replace y, 
by (16) and then integrate over y;, and 2, that 
we have: 


A,(E) (Rk, E)= 


- | 2 
»* ™ 
~ ° 
=~ 
Str 
™ 
— 
7 
6 
~ 
ctr 
6 
~ 
vtr 


N)dl dlodls, (22) 
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where A? = (Le 

The integration ever J» and /; is difficult to 
ko=k;3=0, or, in 
other words, unless the electron is moving in a 


carry out in general unless 
direction perpendicular to the surface. Since the 


exchange potential of an electron depends 
mainly on its velocity, one may obtain a good 
approximation for arbitrary ke, k; by first 
carrying out the integration with ke and k; set 
equal to zero, and then replacing &;, by 
ky?+k2+k,7)! in the final result. If ko=k,=0, 
the integral (22) reduces to: 


xX dre NlE—E v(t), E,)e(h, E)e R,, &). 23 


The integration over \ may be also carried out, 
but before doing this it is best to substitute the 
expressions (21) for g(/;, £) and then integrate 
over £. The integration over \ may then be 
carried out at £=0. If one then divides by the 
value of ¢(k, £) at £=0, one obtains the following 


integral for A,(0). 


e aky jk 
| [3(BC+8l) 
Tp? « km | l 


X log ((Rm?-+/? — kl) / (k —1)?) 
+(/B—kC) tan— (D/(k—)) ] 


+k? log ((B+C+D)/(B+C))}dl (24) 


where B=(p?—k?)!, C=(p?—F)!, D=(k,?—k* 
The last integration may be carried out if p= ~, 
i.e., if the barrier at &=0 is infinite. In this case 


A (0) = (ek, m)[i-1 4a’?+(1/8a*) 
X (1+2a?—3a‘) log ((i+a) 1—a ) |, 


5) 


where a=k;/k,,. It can be shown that this last 
equation is rigorously correct for electrons of 
arbitrary direction if a=(kYP+koe+k3?)! Rm. 

Fig. 2 


potential with @ for p= 


gives the variation of the exchange 
x and ‘for p=1.15k,,. 
The latter were obtained by numerical integra- 
tion. The wave functions for p=1.15k,, corre- 
spond most closely with the computed wave 
functions of the free-electron model. 

These data enable one to get a reasonably good 


picture of the way in which the exchange poten- 
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tials vary near the surface of the metal. The 


type we are considering depend mainly on the 
charge distribution, and not on the exact form 
of the wave functions at the surface. The Fermi 
hole is so large that any local irregularities are 
ironed out when one integrates over the hole to 
obtain the potential. 

The position of the barrier is at &=0 or at 
x=a. The value of a is fixed for any value of p 
by the requirement that the metal be neutral, 
i.e., the requirement of Eq. (4). It may be shown 
by straightforward calculation that this gives: 


(26) 


If p= x, a = 3n, 8k, = 0.61r,. Similarly if 
1.152,,, a=0.11r,. 

If the positions of the barriers are fixed in this 
manner, the exchange potentials are rather in- 
sensitive to the value of p. The values of the 

’ exchange potentials have been computed for 
p=x at a point a distance 0.57, behind the 
barrier, or at x=0.11r,. These values, obtained 

by numerical integration, agreed to within one 
percent with those computed at the same point 
for p=1.15k,,. 

In the actual calculation, the positions of the 
barriers were taken at points slightly different 
from those given above. The wave functions cor- 


responded more closely with those of the free- 


exchange potentials for wave functions of the 
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Fic. 3. Charge distributions for approximate wave fun 


tions and exchange potentials adopted 


with the barriers placed at 


a=0.08r. tor p=«= 


electron model 
a =().58r 


p=1.15k», 


distance 0.03r 


and at and for 
respectively, instead of at points a 
further out. The difference, how- 
ever, is very slight, and hardly enough to affect 
the results appreciably. 

Fig. 3 shows the charge distributions for p=-= 
and p=1.15k, 
for the free-electron model which will be discussed 


as well as the charge distribution 


at more length later on. The charge distributions 
match fairly well, so that it may be expected that 
the exchange potentials will also agree. 

The exchange potentials which have been 
adopted are also shown in Fig. 3. These were 
obtained by drawing smooth curves through the 
known points at x=0.08r, and at 0.587, in such 
a way that they are asymptotic to the known 
values of the exchange potentials in the interior, 
and to zero for large positive values of x. The 
represent the the 
potentials for p=1.15%,, as taken from the curve 
of Fig. 2. The triangles represent the correspond- 


circles values of exchange 


ing values for p= x. One would, of course, expect 
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that the circles are closer to the actual exchange 
potentials than the triangles, but the latter are 
probably not far off. This method for obtaining 
the exchange potentials is, of course, very rough, 
but it was not thought advisable to attempt a 
more accurate calculation at the present time. 
Before discussing the solution obtained, we 
shall, in the next section, discuss the modification 
introduced by the correlation or polarization 


forces which we have so far neglected. 


9. THE CORRELATION ENERGY 


Because of electrostatic repulsion, electrons 
tend to keep as far away from each other as 
possible. In the Slater-Fock approximation, 
electrons of parallel spin keep apart because of 
the exchange effect. There is, however, no allow- 
ance for the interaction of electrons of one spin 
with those of opposite spin except that resulting 
from the integrated charge distribution. The cor- 
relation energy of a metal is the energy gain in 
the binding energy of a metal which results from 
the use of a more general wave function than 
the Slater determinant of wave 
functions. Wigner” computed this energy for free 
electrons by using a wave function which allowed 


one-electron 


for a statistical correlation of the positions of the 
electrons of one spin with those of opposite spin. 
The direct contribution of the correlation hole 
to the work function, computed in (I), is 

. oF) 


ef(r,)—er, f (rs) /3, 27 


where r,f(r,) is the correlation function of 
Wigner.”° 

The correlation 
charge distribution at the surface, and therefore 
the double layer. This latter effect seems to be 
very difficult to take rigorously into account, 


because to do so it would be necessary to abandon 


hole will also influence the 


the one-electron wave functions we have used so 
far. One would expect, however, that the influ- 
ence on the charge distribution could be approx- 
imately taken into account by an appropriate 
modification of the potential curves for the one- 
electron wave functions. This modification would 
be such as to lower the potential in the interior 
of the metal. At large distances from the surface, 


20 E. Wigner, Phys. Rev. 46, 1002 (1934). The correla- 


tion function is represented by the lower curve of Fig. 7. 
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Combined exchange and correlation potentials 
(r,;=4 Bohr units). 


Fic. 4. 


the combined potential of the exchange and cor- 
relation forces should go over into the image 
potential, —e*/4x. 

The amount the potential is lowered in the 
interior will depend on the velocity of the elec- 
tron, and will be smaller for electrons of high 
velocity than for those of lower velocity. In the 
following work we shall make the crude approx- 
imation that the correlation potential varies in 
the same way with velocity and position as the 
exchange potential in the interior and at small 
distances from the surface of the metal, and that 
at large distances from the surface the combined 
correlation and exchange potential approaches 
the image potential asymptotically.” In other 
words, we shall replace A;(x) in Eq. (17) by a 
potential B,(x) where 


B(x) =(1+0)Ax(x) 28) 


except for large x, and we assume that when x 
is large, B,(x) ~e?/4x. We shall take for c the 
ratio of the correlation to the exchange energy 
of the free-electrons. If r,=4 Bohr units, c= 0.24. 

Fig. 4 gives the potentials B;(x) for several 
values of x which were obtained from the ex- 
change potentials, A;(x), of Fig. 3 by the method 
described above. These curves apply for r,=4 
Bohr units.” 


*1 J. C. Slater has suggested a treatment along these lines. 
Cf. reference 8. 

*2 The curves for the smaller values of a should probably 
approach the image potential more rapidly than I have 
indicated. This modification would make little difference 
in the present connection, but should be borne in mind if 
these curves are put to some other use. The author is 
indebted to Professor S!ater for this remark. 
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10. THE SELF-CONSISTENT SOLUTION 


Two solutions of Eqs. (17) have been obtained, 
one using the exchange potentials of Fig. 3, and 
the second using B,(x) of Fig. 4. The first solu- 


tion was very rough, and served only to indicate 





the order of magnitude of the moment of the 
double layer. The Coulomb potential of the 
second was self-consistent to within +0.01e* r,. 
The density of electrons in both cases was taken 
to be approximately equal to the density of free- 
electrons in Na, assuming one free-electron per 
atom, the actual value corresponding to r,=4. 
The solutions were obtained in the following 
manner. One first assumes a charge distribution 
the the 
Coulomb potential, I(x), by integration of Eq. 
2). Eq. (17) is then integrated numerically from 
, large positive values of x into the interior, where 


for electrons, and then calculates 


¢(k, x) must join smoothly with the solution in 
the interior, sin (kjx+-y). This requirement de- 
termines the amplitude of the wave function, and 
also the phase, y, of the plane wave in the 
interior. The integration is carried out for a suf- 
ficient number of values of k,, ke, ky, so that one 
may integrate over the occupied states and 
obtain the charge distribution. If this calculated 
distribution checks the assumed distribution, the 
problem is solved. If not, a new distribution is 
assumed, and the wave functions and charge dis- 
' tribution are redetermined. This process is re- 
peated until a self-consistent solution is obtained. 
The exchange potentials A,(x) (or B,(x)) are 
held fixed throughout. 

In order to determine an accurate self-con- 
sistent solution of Fock’s equations, the exchange 
potentials should also be redetermined from the 
calculated wave functions, but this involves a 
difficult numerical integration and has not been 
carried out. 

In the pure Fock case, the exchange potentials 


¢ A,(x), are used, and the work function, ¢, is 


y= — (h? 2u)Rm? t+el®+Axy,”. (29) 


. Rough calculations that have been carried out 
indicate that the moment of the double layer, 
el”, is about one volt for r,=4. The maximum 
value of the Fermi energy, *k,,?, 2u, is 3.15 ev. 
The exchange energy of the electron of maximum 


velocity is 0.61 e r,=4.1 ev. Thus if we neglect 
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WORK FUNCTION 
the correlation energy entirely, the work function 
will be about two volts. If we include the direct 
contribution of the 


section 8) the work function is about three volts. 


correlation energy (cf 

It is not altogether reasonable to include the 
direct contribution of the correlation energy if 
we omit it in the calculation of the surface double 
layer (V°) because the work function is rather 
insensitive to the exchange potential employed. 
Our neglect of the correlation energy has been 
we had 
still 


partly compensated by a large 1°. If 
omitted the exchange energy, we would 
obtain a positive (though probably small) value 
for the work function with a very large V° (~4 
volts). 

The solution obtained when the potentials 
B,(x) of Fig. 4 were used is illustrated in Fig. 5. 
The potential curves shown include the Coulomb 
potential, I(x). The moment of the double layer, 


170 


1, is 0.06e? r,, or about 0.4 volt. The value of 
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Fic. 5. Self-consistent charge®distribution and potentials 
for the free-electron model (r,=4 Bohr units). The po- 
tential curves shown include the exchange and correlation 
potentials of Fig. 4 and the Coulomb potential of the 
charge distribution. 
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the work function is then 2.35 volts, if we include 
the direct contribution of the correlation energy 


(1.0 ey 


11. COMPARISON WITH SODIUM 


It has been shown by Wigner and Seitz? and 
by Slater® that the wave functions of the con- 
ducting electrons of the sodium crystal are, except 
in the immediate vicinity of the ions, very nearly 
the same as those of perfectly free electrons. The 
density of electrons is practically uniform 
throughout the crystal. The Fermi energy is 
within one or two percent of the classical free- 
electron value, 3h°k,,2, 10u. The exchange poten- 
tials are very nearly constant in space, and are 
close to the free-electron values in magnitude. 
One would thus expect that the free-electron 
model which has been discussed in the preceding 
sections would bear a close relationship with the 
sodium crystal. 

Since the exchange potentials are constant in 
the interior of the lattice, the Fock equations are 
identical with the Hartree equations. The ex- 
change terms affect the energies of the electrons, 
but not the wave functions. If we omit the 
exchange terms, the potential field in which the 
electrons move is the combined potential of the 
ions and the Coulomb potential of the electrons. 
Wigner and Seitz replace the actual (body- 
centered) lattice by a set of close packed spheres. 
They found that if the zero of potential is taken 
on the surface of an s sphere, the energy of the 
lowest state is Ej = —0.23 ev. The actual poten- 
tial on the surface of an s sphere, relative to a 
point outside one of the faces of the crystal, wil 
differ from zero by the moment of the double 
layer of this face. 

The work function, g, is (cf. (1): 

Ona = — Lo teDt+A in? +C—hPk,2/2u, (30) 
where we have written C for the contribution 
from the correlation energy, and eD for the 
moment of the double layer. The only difference 
between (30) and the corresponding formula for 
the free-electron model is that — E)+eD replaces 
el” of the latter. 

One might at first be tempted to identify e|”° 
with eD. Some care must be taken, however, as 
the definitions of the moments in the two cases 


TABLE I. Comparison of theoretica lues of the 
function of the free elk [ n model wit perimental valu 
for 

FR FRE 
ELECTRO I 
MODEL, Mop 
CORRELA- CORRE A > 
1ON FORCES TION Fo DI 
OmMITTI Inct ) (ex 
Moment of double laver 1.0 0.4 0.15 
Work function 2.0 2.35 2.46 
> 26 
* A. R. Olpin, see reference 11, p. 75 


t Z. Berkes, Math. Phys. Lapok 41, 131 (1934 


are slightly different. Thus eD will be zero if the 
distribution of electrons in the surface s spheres 
is the same as that in the inner s spheres, while 
e\° will be zero if the distribution is constant for 
x <0 and zero for x >0. For the same distribution 
of electrons, eD will be somewhat smaller than 
el. The difference can be computed for any face 
of the crystal, and it amounts to approximately 
0.25 ev for the 110-face. Values for other faces 
are slightly larger. 

Thus —Ey+eD for the 110-face should be 
about the same as eJ”, while the corresponding 
values for other faces should be a little less. 
This, of course, rests on the assumption that the 
electronic charge distribution at the surface of 
the sodium crystal is approximately equal to 
that of the free-electron model. Calculations that 
have been attempted for the 110-face indicate 
that this assumption is justified. 

Table I summarizes the results that have been 
obtained for the free-electron model, together 
with the rather uncertain experimental values of 
the work function of sodium. The agreement is 
very satisfactory, considering the approximate 
nature of the theory. The moment listed for 
sodium was obtained by reducing the value for 
the free-electron model (0.4 volt) by 0.25 volt. 


12. CONCLUSIONS 


The following are the main results obtained 
from the present investigation of the free-electron 
model : 

(1) The barrier at the surface is due largely to exchange 
and polarization forces, rather than to ordinary electrostatic 
forces. 

(2) There is no single potential barrier which is satis 
factory for electrons of all velocities. Eaci: electron has, s 


to speak, its own barrier. 
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(3) The moment of the double layer at the surface is 
relatively small, in agreement with the conclusions reached 
in reference I. This, of course, means that the ordinary 
electrostatic potential in the interior is small. 

(4) Approximate agreement is obtained between the 
work function of the free-electron model and the experi- 
mental values of the work function of sodium. 


The writer wishes to express his gratitude to 
Professor E. Wigner, under whose direction this 
work was carried out, and to Professor J.. H. 
Van Vleck, Professor J. C. Slater, and Dr. F. 
Seitz for interesting discussions. 


APPENDIX 


Evaluation of the integral: 
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dy dz}. 


7! (E&P +0 —-7nF Tes, 
This integral is not absolutely convergent, and in evaluat- 
ing it one must be careful to introduce polar coordinates, @, 
p in place of y; and 2, and then integrate first over 6. If 
one chooses the origin of @ properly, the integral takes the 
form: 


ha he - 06 


eo . eidp 
I= lp dé 
J ol Gar 


with \2 = (ke—/2)?+(ks—/;)*. The integration over @ gives 
Bessel’s function Jo(Ap). We thus have: 


wes) 
I=2 | 


Eq. (22) follows immediately. 


pJo p)dp 
((E—£1)?+p*)! 
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The Variation of Young’s Modulus with Magnetization and Temperature in Nickel 


SIDNEY SIEGEL AND S. L. QuimBy, Columbia University 
(Received March 3, 1936) 


The relation between Young’s modulus, magnetization, 
and temperature in annealed polycrystalline nickel, 99.7 
percent pure, is determined for values of J between zero 
and saturation, and for values of T between 23.5°C and 
400°C. The percent increase in Young’s modulus is propor- 
tional to J? between zero and about 0.4 saturation, at all 
temperatures below 311°C. The total increase in the 
modulus from the demagnetized to the saturated state is 
6.7 percent at 23°, reaches a maximum of 18.7 percent at 
185°, and decreases to zero at the Curie point. The results 
indicate that the theory of this phenomenon offered by 


HE domain theory of ferromagnetism, first 

suggested by Webster in connection with 
magnetostriction and since developed by Heisen- 
berg, Akulov and others, is noteworthy for its 
success in correlating different ferromagnetic 
phenomena in both single and polycrystalline 
materials. For example, Akulov' has shown that 
the change of Young’s modulus with magnetiza- 
tion in a polycrystalline substance can be de- 
scribed solely in terms of the initial susceptibility 
and one magnetostrictive constant of a single 
crystal. The present paper is a report of an 
experimental investigation of the relation be- 
tween Young’s modulus and magnetization in 
annealed polycrystalline nickel at various tem- 





' Akulov, Zeits. f. Physik 85, 661 (1933). 


Akulov is probably essentially correct but requires modi- 
fication. The method employed for the measurement of 
Young's modulus yields at the same time a measure of the 
coefficient of internal friction of the material. The internal 
friction decreases with increasing magnetization at all 
temperatures below the Curie point. In the demagnetized 
material it reaches a maximum value at the same temper- 
ature as does the change in elastic modulus. Its value at 
magnetic saturation is the same order of magnitude as that 
of a nonferromagnetic substance. 


peratures below the Curie point. The results are 
compared with the formulae of Akulov’s theory. 


EXPERIMENTAL METHOD 


The research consists in ascertaining the value 
of the Young’s modulus, £, corresponding to an 
intensity of magnetization, J, and a temperature, 
T, for values of J lying between zero and tech- 
nical saturation, and values of T ranging from 
20°C to 400°C. Accordingly the experimental 
method must permit the simultaneous measure- 
ment of EZ, J, and T. To this end a preliminary 
determination is made of the relation between 
J, T and the normal magnetic induction, B, 
so that in the subsequent procedure a knowledge 
of B and T suffices to fix the value of J. 
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Fic. 1. Cross section of the assembly. 


A specimen is prepared in the form of a right 
circular cylinder 0.464 cm in diameter and 7.2 cm 
long terminated by the end sections of a coaxial 
ellipsoid of revolution. The form approximates 
that of an ellipsoid of revolution whose major 
axis is 17.8 cm and minor axis is 0.464 cm. The 
specimen is mounted axially in a tubular electric 
furnace surrounded by a solenoid, and the 
variation of the normal magnetic induction with 
temperature and applied magnetizing field, //o, 
measured with a fluxmeter galvanometer. 
Corresponding values of J are then calculated 
with the formula J = (B—J/))/(4%—N), where NV 
is the demagnetization factor of the ellipsoid, 
assumed the same as that of the specimen. 
The latter assumption is justified by the ob- 
served uniformity of B over the length of the 
cylinder, and by the fact that, since V=0.028, 
an error in the first significant figure would alter 
the calculated value of J by less than 0.1 percent. 
After these measurements are completed the 
cylindrical part of the specimen is sawed out 


1S 


and its ends ground flat. 

The method employed for the measurement of 
Young’s modulus is a modification of that 
previously described in this journal by Zacharias.’ 
The specimen cylinder of nickel forms one part 
of a composite piezoelectric oscillator con- 
structed by cementing to one end a cylinder of 
crystalline quartz of square cross section equal 
in area to that of the specimen. The quartz 
cylinder is so cut from the crystal that the optic 
axis lies perpendicular to one pair of opposite 
faces and an electric axis to the other, and the 
latter pair is coated with gold leaf. An har- 
monically varying potential difference of con- 





2 J. Zacharias, Phys. Rev. 44, 116 (1933). Other papers 
on the use of composite piezoelectric oscillators for the 
measurement of elastic constants are: L. Balamuth, Phys. 
Rev. 45, 715 (1934); and F. C. Rose, Phys. Rev. 49, 50 
(1936). 
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Fig. 2. Distribution of magnetic flux in the specimen. 
The dotted curve shows the distribution when the speci- 
men is placed in a long solenoid. 


stant amplitude is established between the gold 
leaf electrodes, and means are provided for 
measuring the amplitude of the resulting current. 
In consequence of the harmonically varying 
piezoelectric stress in the quartz, produced by 
the electric field between the electrodes, the 
entire oscillator assumes a state of forced longi- 
tudinal vibration, the frequency of which is the 
same as that of the applied potential difference. 

The amplitude of the current which flows to 
the quartz varies critically with frequency in the 
neighborhood of certain ‘‘resonance frequencies”’ 
at which the amplitude of the vibration passes 
through a maximum. If fo denote one of these 
frequencies, then at a frequency f’ slightly less 
than fy the current passes through a maximum, 
and at a frequency f” slightly greater than fy it 
passes through a minimum. It shown in 
Zacharias’s paper that if J’ and J” 
respectively, the maximum and minimum values 
of the current amplitude, then fy is given by the 
formula* 


fo=f' +l” (f"—-f') (7’+]"’). 


Furthermore, the fundamental frequency, f;, of 
free longitudinal vibration of the nickel cylinder 
alone is related to fy by the formula 


mif, tan (afo/fi)+mofe tan (rfo/ fe) =0, 


where m, and my, are the masses, respectively, of 
the nickel and quartz cylinders, and f2 is the 
fundamental frequency of free longitudinal vibra- 
tion of the quartz cylinder alone. The value of fs 


3 Misprints which occur in Eqs. (11) and (12) of Zach- 
arias’ paper are here corrected. 


is 


denote, 


(1) 


(2) 














is obtained by observing the current which 
flows to the quartz without the specimen at- 
tached. Lastly, the quantity f; is related to the 
Young’s modulus of the material by the formula 





E=4Pf 2p, (3) 


where / is the length of the specimen cylinder and 
p is its density. 

This method for measuring Young’s modulus 
yields at the same time the value of the loga- 


rithmic decrement, 6, of free longitudinal vibra- 
tion of the specimen cylinder. Thus 


6 = 2r(1+me2/m,){(f" —fo)(fo—f’) | */fo. (4) 





The coefficient of internal friction of the material 
is equal to 6, E. 
Errors introduced by the adhesive and by 
’ inequality of areas at the interface are minimized 
by adjusting the lengths of the cylinders so that 
fi=fe within 10 percent. These frequencies lie in 
the neighborhood of 36 kilocycles. 

Fig. 1 is a cross-sectional diagram of the 
assembly upon which the measurements are 
made. The composite oscillator is mounted in a 
thin glass tube, G, upon an annular duralumin 
knife-edge, K,, and a straight knife-edge, Ko, 
located at displacement nodes of vibration. The 
temperature of the specimen is measured with a 
Chromel-P vs. Alumel thermocouple in contact 
) with the knife-edge K,. The magnetic induction 
is measured with a fluxmeter galvanometer con- 
nected to the coil C;, constructed by winding 90 
turns of No. 40 chromoxide covered copper wire 
on a thin cylinder of baked lavite. The mag- 
netizing coils Co, C3, and Cy, are designed to 
secure as nearly as possible a uniform mag- 
netization over the length of the specimen 
cylinder. C2, and C; are each 92 turns of copper 
ribbon 0.023 cm thick and 0.635 cm wide, in- 
sulated with thin strips of mica. C, is 340 turns 
of No. 20 chromoxide covered copper wire wound 
in 6 layers 5 cm long on a Pyrex tube 1 cm in 
diameter. The layers of this coil are also sepa- 
. rated by mica sheet. The magnitudes and direc- 

tions of the end and middle coil currents are so 


adjusted that B is uniform to 5 percent over 
eight-tenths the length of the specimen, at all 
temperatures and for all values of B. Typical 


distribution curves for B are shown in Fig. 2. 
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The assembly shown in Fig. 1 is placed in a 
rectangular duralumin oven, the interior dimen- 
sions of which are 7 in. X7 in. X12 in. The oven 
is heated electrically over its entire outer surface 
with grids formed by laying No. 22 constantan 
wire along parallel grooves cut in sheet transite. 
The outer thermal lagging is corrugated asbestos 
paper built up to a thickness of 1.75 in. 

The voltage applied to the quartz crystal is 
supplied by a vacuum tube oscillator whose 
frequency can be varied in the neighborhood of 


fo by means of a variable condenser. As noted 


above, the determination of E requires the 
measurement of the maximum and minimum 
piezoelectric oscillator current amplitudes, and 
the frequencies at which these values occur. 
In previous applications of the method, the 
recording device for measuring the currents has 
been a d’Arsonval galvanometer, used in con- 
nection with a vacuum thermocouple or other 
rectifier. In the present instance, however, the 
heat produced by the current in the magnetizing 
coils changes the temperature of the specimen 
too rapidly to permit the use of this instrument. 
The entire set of observations must be made 
within six seconds after the magnetizing field is 
applied. This is accomplished in the following 
manner : 

The difference of potential across a non- 
reactive resistance in series with the crystal is 
amplified and rectified with a vacuum tube 
rectifier. The rectified voltage is, in turn, ampli- 


Fic. 3. Polar oscillogram of the piezoelectric oscillator 
current amplitude. 





666 S. SIEGEL AND 
t | Nie we L5G 
ou QUARTZ ” b 
as r oS 0 rH \ 
20 S 
<> 3 6 ' 
| 3 . +25 +1385 ~-9 +90 «+135 





if 








=> 
> 








8s> 
0 ‘ 
0 ° io § 
t = 
“ = 
+135 
NW 0.01 0.001 
= mh libs 
45S 4.5 


Fic. 4. The amplifying circuit. Resistances are in megohms 
and capacities in millimicrofarads. 


fied and energizes a General Electric Company 
Type PM-11 supersensitive oscillograph element. 
The light beam from the oscillograph passes 
through a cylindrical lens and falls upon a sheet 
of bromide paper attached to a vertical circular 
disk. This disk is mounted centrally on one end 
of the horizontal shaft of the variable condenser 
which controls the exciting frequency. An electric 
motor rotates the condenser at a rate which 
carries the frequency through f’, fo, and f” in 
4 seconds. During this time the light beam traces 
on the sensitized paper a polar diagram of the 
piezoelectric current amplitude, in which the 
angular coordinate is proportional to frequency 
and the radial coordinate to current. The con- 
denser is frequency calibrated by comparison 
with the output of a piezoelectric clock,‘ and 
frequency calibration marks are placed on the 
record by illuminating it locally for an instant 
through short diametrically opposed radial slits.*® 
A typical oscillogram is reproduced in Fig. 3. 
Fig. 4 is a diagram of the amplifying circuits 
employed. The rectified voltage corresponding to 
an oscillator current amplitude between J’ and J” 





‘Quimby, Phys. Rev. 39, 345 (1932). 

5 It should, perhaps, be noted that this polar oscillograph 
is made fully automatic in operation by employing the 
various features of an RCA model 86 remote control 


radio unit. The different station selector contacts serve to 
start and stop the motor, flash the light source, open the 
source slit, etc. 
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Fic. 5. Curves showing the resonance frequency (circles) 
and decrement (dots) of the piezoelectric oscillator as 
functions of magnetization at 177°C. 


is balanced out with the 10,000-ohm _ potenti- 
ometer, so that only the fluctuation associated 
with J’ and J” is recorded by the oscillograph. 
The deflection of the oscillograph mirror varies 
linearly with the amplitude of the potential 
difference across the resistance, R, in series with 
the oscillator, and so with that of the current. 
Furthermore, the current is independent of the 
magnitude of R, since the latter is small in com- 
parison with the impedance of the oscillator. 
Accordingly each oscillogram is current cali- 
brated by noting the values of R required to 
duplicate the deflections corresponding, respec- 
tively, to I’ and J”, the current being held con- 
stant at the intermediate value. The actual 
working formula for the calculation of fo is thus, 
from Eq. (1), 


fo=f' +R" (f"—f’) (R’+ R”), 


where R’ and R” correspond to the deflections 
associated, respectively, with J’ and J”. 

The experimental procedure is as follows: The 
temperature of the oven is stabilized, the speci- 
men is placed in the desired state of magnetiza- 
tion on the normal induction curve by repeated 
reversals of the magnetizing current, and the 
resonance current record is taken immediately. 
The flux density, B, is then measured, the oven 
is allowed to resume its former equilibrium 


(5) 
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temperature, and the process is repeated for 
various values of J. The oscillograph records are 
reduced with the aid of Eqs. (5) and (4), and 
curves are plotted, of which those shown in 
Fig. 5 are typical, showing fy) and 6 as functions 
of J/J,,, where J,, denotes the saturation mag- 
netization at the temperature of observation. 
Lastly, isothermals of this sort are obtained at 
various temperatures. 


RESULTS 


The specimen material upon which the obser- 
vations here reported were made is Mond nickel 
containing 0.295 percent impurity distributed as 
follows :® Co, 0.00; Si, 0.01; Cu, 0.02; Fe, 0.11; 
Mn, 0.00; C, 0.049; Mg, 0.11; S, 0.005. The 
specimen was annealed at 1100°C for 4 hours in 
hydrogen at atmospheric pressure and allowed 
to cool to room temperature in 18 hours. Micro- 
scopic examination reveals it to be an array of 
crystals whose average linear dimensions are 
0.4 mm. The coercive force is 0.3 oersted, and the 
remanent intensity of magnetization is in the 
neighborhood of 9 gauss. The Curie point is 
353°C, as determined to within one degree by the 
magnetic measurements and by the change of 
Young’s modulus with magnetization and tem- 
perature. The variation of /,, with temperature 
is given by the first two columns of Table I. 

The variation of Young’s modulus with mag- 
netization and temperature is shown by the 
curves of Fig. 6. The points on these curves were 
calculated with Eqs. (2) and (3), and values of fo 
secured from smooth curves of the sort shown in 
Fig. 5. 

Table II gives the variation with magnetiza- 
tion and temperature of the quantity AE/E 
defined by the formula, AE/E=(E ,—Eo)/Eo. 


TABLE I. Variation of J,, with temperature. 





1 te Eo X10 7 [* Eo X10"! 
23.5 485 20.82 270 332 18.54 
70.0 472 20.09 311 260 19.10 

106 459 19.26 340 162 19.63 
147 440 18.37 353 39 19.79 
177 422 17.93 365 2 19.70 
197 408 17.81 369 1 19.68 
200 405 17.78 385 - 19.57 
217 371 17.85 397 


19.49 








_ *This analysis, together with the specimen material 
itself, was furnished by the International Nickel Company. 
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E, and £» denote the values of Young’s modulus 
when the value of the magnetization is, respec- 
tively, J and zero. Ep is given as a function of 
the temperature in the third column of Table I. 

The variation of the logarithmic decrement 
with magnetization and temperature is shown by 
the curves of Fig. 7. 


ACCURACY OF THE MEASUREMENTS 


Temperature 

The Chromel-P vs. Alumel thermocouple was 
calibrated before and after the measurements 
against the melting points of ice, Sn, Pb and Zn. 
0.25°C. The 
temperature variation in the oven over the region 


The two calibrations agreed to 
occupied by the specimen is less than 0.3°C. 


That along the specimen is, of course, con- 
siderably less. The temperature of the thermo- 
couple in contact with the knife-edge (Fig. 1) 
agrees with that of the specimen to better than 
0.1°C. Accordingly the error in the temperature 


measurement is not greater than 0.5°C. 


Magnetic 

The measurement (on the approximate ellip- 
soid) of the saturation magnetization is accurate 
to 1 Values of J between 
saturation are reliable to 2 percent, and errors 


percent. zero and 
arising from this uncertainty are reduced by the 
smoothing process represented by the curves of 
Fig. 5. 
Elastic 

The relation and fi 
involves the length and density of the specimen 
cylinder at the observation. 
These are evaluated with accuracies of 0.02 and 
0.05 percent, respectively, from measurements 


(3)) between E 


(Eq. 


temperature of 


TABLE II. (Ey—Eo)/Eo. 
Wd 

> 0.2 0.4 0.6 0.8 1.0 
23.5 | 0.0040 0.015 0.035 0.055 0.067 
70.0 .0065 025 052 .077 .096 
106 0088 .033 071 .109 131 
147 .0109 .048 .094 .136 .170 
177 .0084 .046 .106 .156 .186 
197 013 .053 117 .163 .186 
200 .013 055 .116 .161 .188 
217 .013 .053 111 .152 .176 
270 .0092 .039 .073 .098 .109 
311 .0063 .023 .046 .055 .058 
340 .0025 .0061 O11 013 .016 
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Fic. 6. Curves showing the variation of Young’s modulus 
with magnetization and temperature in nickel. The num- 
bers on the curves give the magnetization in fractions of 
the saturation value. 


made at 30°C, together with values of the 
coefficient of thermal expansion of nickel re- 
ported by Williams.’ 

A possible source of error in the measurement 
of Young’s modulus lies in the fact that about ten 
percent of the length of the specimen cylinder at 
each end is in a lower state of magnetization 
than that observed at its center. The maximum 
effect of the consequent elastic inhomogeneity 
on f,; is readily obtained by an analysis similar to 
that which leads to Eq. (2). It follows that if 
the value of E over the end sections is assumed 
to be 0.9 the value over the central eight-tenths 
of the cylinder, then the resulting error in the 
calculated value of E is only 0.25 percent. This 
is merely a consequence of the fact that the 
fundamental frequency of free longitudinal vibra- 
tion of a rod is principally determined by the 
inertia of its terminal parts and the elasticity of 
its If the above-mentioned in- 
homogeneity occurred at the of the 
cylinder the resulting error in E would be 4.2 


central part. 


center 


percent. 

The major part of the uncertainty in E arises 
from the flatness of the maxima and minima 
shown on the oscillograms when the decrement, 
6, is large. The error thus introduced into the 
determination of f’ and f” can affect the calcu- 
lated value of E by as much as 0.2 percent. 


7 Williams, Phys. Rev. 46, 1011 (1934). 
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Fic. 7. Curves showing the variation of logarithmic 
decrement with magnetization and temperature in nickel 
The numbers on the curves give the magnetization in 
fractions of the saturation value. 


Finally, the reproducibility of the results is 
indicated by the observations at 197° and 200° as 
plotted in Fig. 6. The second set were taken a 
month later than the first. As a further test, the 
oscillator was demounted, disassembled, re- 
assembled, and remounted, and a new set of 
data was taken with the specimen demagnetized 
and magnetized to saturation. The lowest and 
highest curves of Fig. 6 were reproduced with 
maximum accidental errors in E of 0.3 percent 
and 0.16 percent, respectively. 

The values of 6 obtained by the present 
method are reliable to 10 percent. 


THEORETICAL DISCUSSION 


The domain theory of ferromagnetism is based 
on the assumption that a ferromagnetic crystal 
is composed of a large number of regions, each 
of which is spontaneously magnetized to an 
intensity which is a function of the temperature 
alone. When the holomagnetization is zero these 
micromagnetizations are aligned in equal num- 
bers along the various directions of ‘‘easy mag- 
netization” of The effect of a 
magnetizing field is, in the first instance, simply 
to alter this uniform distribution, that a 
preponderant number of micromagnetizations lie 
in that direction of easy magnetization which 
most nearly coincides with the field. When this 


the crystal. 


so 


process is complete a further increase in the field 
strength of the micro- 
magnetization of each domain into the direction 
of the field. 

Each domain is deformed in a manner which 


produces a rotation 

















YOUNG’S MODULUS AND M 
depends upon the direction of its magnetization. 
Consequently any alteration in the directional 
distribution of the micromagnetizations through- 
out the crystal is accompanied by a deformation 
of the crystal as a whole. When the alteration is 
produced, as above, by a magnetizing field alone, 
the resultant deformation constitutes the ordi- 
nary phenomenon of magnetostriction. 

Now it is well known that the holomagnetiza- 
tion of a ferromagnetic body in a given mag- 
netizing field is, the 
application of a stress, and furthermore, that 
the magnitude of this alteration depends upon the 
magnetization, being nil when the substance is 


in general, altered by 


magnetized to saturation. In terms of the domain 
theory this means that a given applied stress 
alters the distribution of the micromagnetiza- 
tions, and hence the linear dimensions of the 
body, by an amount which depends upon the 
holomagnetization. This is the qualitative expla- 
nation of the phenomenon dealt with in this 
paper. 

The problem is attacked quantitatively by 
Akulov,' who employs the analysis invented by 
Heisenberg* to resolve the simpler problem of 
single crystal magnetostriction. The procedure is 
simply to calculate the most probable distribu- 
tion of the micromagnetizations for a given 
applied field and stress, subject to the conditions 
that the total number of and the 
quantity U are constant, where U is defined by 


domains 


the formula 


U= — (11,J,+HeJ2+H3J3) — Fx. 


In this expression 7; and J; are the components 
of the field and magnetization respectively, F is 
the applied force per unit area, and X is that 
part of the change in length per unit length 
parallel to F which depends upon the distribution 
of micromagnetizations. The calculation yields 
the distribution in terms of 7, J, F, and the 
magnetostrictive coefficient of a domain for a 
direction of easy magnetization. \ is evaluated in 
terms of the same quantities by summation over 
the array of domains, and from this the variation 
of the elastic modulus with magnetization. The 
behavior of a polycrystalline material is obtained 
by averaging over all directions in the crystal. 
The calculation cannot be completed without 





* Heisenberg, Zeits. f. Physik 69, 287 (1931). 
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Fic. 8. AE/Epo is proportional to the square of the mag 
netization for values of J up to 0.4 saturation. 


approximation. Akulov obtains two formulae, 
one of which is valid when the applied field is 
small, and the other when the field is sufficiently 
large that the effect of stress on the distribution 
is nil. Thus, for a polycrystalline specimen of 


nickel in which the distribution of micro- 
crystalline axes is random, 
E;-—Eo 81 xo 
- Aur(J/J..)*, (0) 
Ey 35 J. 


provided J/J,, is small, while if J=/J,,, 
E. — Ey Dx oA Eo i5J,.? 


Ey 1—9yoAu11%Eo/15J,2 


where E,, denotes the value of Young’s modulus 
when the material is magnetically saturated. 
In arriving at these expressions /7 has been ex- 
pressed in terms of J by the formula, J = xo/7. 
Hence xpo is the initial magnetic susceptibility of a 
single crystal of nickel. A,,; is the saturation 
magnetostrictive coefficient for the direction of 
easy magnetization in the crystal. The measure- 
ments of Masiyama® give Ai1=27X10~° and 
those of Becker and Kersten'® give xo= 20, both 
at room temperature. These formulae and data 
will now be compared with the results of the 
present research. 


® Masiyama, Tohoku Imp. Univ. Sci. Rep. 17, 945 
(1928). 
‘0 Becker and Kersten, Zeits. f. Physik 64, 660 (1930). 





670 J. J. BRADY AND 

The straight lines shown in Fig. 8 indicate 
that, in accordance with Eq. (6), (E,y—Eo)/Eo 
varies linearly with (//J,,)? up to temperatures 
at least as high as 311°C, and for values of J/J,, 
as large as 0.4. 

The calculated value of the of 
(J/J,,)? in Eq. (6) is 0.30. The observed value at 
23.5°C is 0.10. The value of (E,,—Eo)/Eo calcu- 
lated from Eq. (7) is 0.085. The observed value 


coefficient 


at the same temperature is 0.067. An approxi- 
mate test of the correctness of Eqs. (6) and (7), 
independent of the values of xo, A111: and E can be 
obtained by noting that the value of the nu- 
merator on the right-hand side of Eq. (7) is only 
0.078. Hence the ratio of this quantity to the 
coefficient of (J/J,.)? in Eq. (6) should equal 
35/135, or 0.26, with an accuracy of 8.5 percent 
whatever the value of xo, Ain, and FE. The ob- 
served ratio is 0.67. 

These discrepancies cannot be ascribed to a 
nonrandom distribution of microcrystals in the 
specimen. The possible values of the ratio noted 
in the preceding paragraph range from 0.20 for a 
(111) direction in a single crystal to 0.33 for 
a (110) direction. Nor can they be ascribed to 
the neglected effect of temperature, for it is 
observed that the ratio decreases with increasing 


temperature from the value 0.67 at 23.5°C to 
the value 0.39 at 311°C. 
MAY 1, 1936 PHYSICA 
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von Auwers"! has investigated the AE effect in 
nickel at room temperature by a method dif- 
ferent from that here employed. He obtains the 
values 0.30 and 0.43 for the ratio, as observed on 
specimens characterized, respectively, as “‘soft’’ 
and “‘very soft.”’ 

In conclusion, therefore, it is evident that, 
while the domain theory is probably essentially 
correct, the greatly simplified calculation of 
Akulov requires emendation of such nature as to 
alter at least the numerical coefficients in Eqs. 
(6) and (7). 

The interpretation of the decrement measure- 
ments awaits the development of a theory of 
internal friction in ferromagnetic substances. It 
may be noted that the observed decrement of the 
magnetically saturated nickel (curve e of Fig. 7) 
is the same order of magnitude as that of non- 
ferromagnetic materials. 

The authors gratefully acknowledge their in- 
debtedness to Mr. C. A. Crawford of the Inter- 
national Nickel Company for the specimen 
material and its analysis; to Mr. F. A. Nealon 
of the American Enameled Magnet Wire Com- 
pany for the gift of the chromoxide covered wire ; 
and to Professor William Campbell and Mr. 
Maurice Pinel of Columbia University for their 
assistance in the microscopic study of the poly- 
crystalline specimen. 


yon Auwers, Ann. d. Physik 17, 83 (1933). 
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Photoelectric Properties of Sodium Films on Aluminum 
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Previous work has shown that thin films of potassium, 
rubidium, and caesium deposited on a silver surface give a 
maximum sensitivity when the thickness of the film is of the 
order of ten molecular layers or less. As the deposit con- 
tinues the sensitivity decreases before reaching a constant 
value. Sodium films were deposited on aluminum, and it 
was found that the maximum sensitivity was not attained 
until the thickness of the film had reached a value of 
eighty molecular layers. The sensitivity was then found to 
remain constant for greater thicknesses. No detectable 
photo-current was observed until the film reached a thick- 
ness of five molecular layers. Spectral distribution curves 
were obtained and analyzed according to Fowler's theory. 
This yielded a value of 5150A for the threshold of the 


sodium film. These observations indicated that the thresh- 
old was the same for all thicknesses of film. Current-voltage 
curves were obtained for different thicknesses of film, and 
these were analyzed according to DuBridge’s theory. These 
results also indicate that the work function of the various 
films is the same. It was concluded from these data that the 
work function of the freshly distilled aluminum on which 
the deposits were made was 4.08 volts. This gives a value of 
3020A for the threshold of aluminum. It was found that the 
photoelectric current decreased immediately after forming 
a film, but that the decrease was much less for films greater 
than fifty molecular layers. It is suggested that this de- 
crease results from the surface migration of the sodium 
atoms in forming aggregates. 
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INTRODUCTION 


T is well known that when only a small fraction 
of the surface of a tungsten or platinum plate 

is covered with either Cs, Rb, or K atoms, the 
work function of the composite surface is much 
less than that of the clean underlying metal. 
Langmuir and Kingdon,! and also Becker,” have 
taken the view that atoms may be adsorbed as 
ions on surfaces if the work function of the 
surface is greater than the ionization potential 
of the adsorbed atom. According to Becker's 
“ion-grid”’ theory, ions distributed here and there 
over the surface of a metal lower the work 
function of the composite surface. The criterion 
stated above, for the existence of adsorbed ions, 
must be taken only as a rough guide. Gurney*® 
has pointed out, from quantum-mechanical con- 
siderations, that the criterion could be expected 
to apply only to single valence atoms (e.g., 
alkali metals) condensed on a surface. Experi- 
mental results for deposits of alkaline earth 
metals are in agreement with these views.’ In 
terms of the picture discussed by Gurney, the 
criterion for the existence of an ion grid on the 
surface should be stated in terms of probability. 
If the work function of the underlying surface is 
greater than the ionization potential of the 
adsorbed atom, the probability of its being ad- 
sorbed on the surface as an ion is quite high. 
On the basis of this theory, there is a finite 
probability for the existence of adsorbed ions on 
the surface even though the ionization potential 
of the deposited atoms is equal to or even 
slightly higher than the work function of the 
underlying metal. The probability would be 
expected to decrease to a very small value if the 
ionization potential of the atom were much 
greater than the work function of the surface. 

It was the purpose of this investigation to 
choose a combination such that the ionization 
potential of the deposited atoms was considerably 
higher than the work function of the underlying 
metal. It was desirable also to choose a combina- 
tion that would give a photoelectric effect with 


‘1. Langmuir and K. H. Kingdon, Proc. Roy. Soc. A107, 
61 (1925). , 

2]. A. Becker, Phys. Rev. 28, 341 (1926). 

*R. W. Gurney, Phys. Rev. 47, 479 (1935). 

‘J. A. Becker, Trans. Faraday Soc. 28, 155 (1932). 
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visible light. Sodium films deposited on alu- 
minum satisfied these conditions. 

Brady® reported in studies of the photoelectric 
properties of K, Rb, and Cs films deposited on 
silver that the maximum sensitivity was reached 
when the calculated thickness of the film was in 
the order of ten molecular layers or less (depend- 
ing on the frequency of the exciting light). 
In the same report a curve was shown for sodium 
deposited on tungsten and the maximum sensi- 
tivity was not reached until the thickness of the 
film was well over one hundred molecular layers. 

Fig. 1 is an unpublished curve from Brady's 
data and it is interesting to contrast this with 
Fig. 2. 


EXPERIMENTAL PROCEDURE 


The photoelectric cell was almost identical to 
that used in a previously reported experiment.® 
Likewise, the same procedure was used for 
evacuating the cell. After a thorough outgassing 
treatment, the cell was sealed off from the pumps. 
During the course of an experiment, the charcoal 
trap attached to the cell was immersed in liquid 
air. This reduced the pressure within the cell to 
less than 10-§ mm of Hg as registered by the 
ionization gauge. The sodium was outgassed by 
repeated distillation before being driven into the 
lower part of the photoelectric cell. The alu- 
minum which formed the underlying surface was 
distilled from 30 mil tungsten spirals. 

All photoelectric studies were made immedi- 
ately after depositing a fresh coat of aluminum 
on the inside surface of the cell. This was done 
to insure an outgassed surface on which to 
deposit the sodium films. 

A mercury arc with suitable filters was used 
as a light source in most of the experiments. 
The spectral distribution curves were obtained 
by using a tungsten filament lamp in conjunction 
with a double prism monochromator of the van 
Cittert type. 

As in the previous work,’ a monomolecular 
layer is defined on the assumption that the 
atomic spacing for the sodium atoms is the same 
in the film as in a solid mass of the metal. It is 
assumed, also, that the underlying surface is 
perfectly smooth. On these assumptions, a mono- 


° J. J. Brady, Phys. Rev. 41, 613 (1932). 


6 J. J. Brady, Phys. Rev. 46, 768 (1934). 
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Fic. 1. Variation of photo-current with the thickness of a 


potassium film deposited on silver (after Brady). 


molecular layer of sodium is equivalent to a 
deposit of 8.7 X10" atoms/cm’. 


RESULTS 


Fig. 2 shows a typical curve for the variation 
of photo-current with the thickness of the film. 
No photo-current was detected before the dep- 
osition of the sodium, even in the presence of a 
very intense white light. The Pyrex window of 
the cell was effective in cutting out the shorter 
wave-lengths which would cause the aluminum 
to emit photoelectrons. The photo-current from 
the sodium film became detectable only when the 
thickness of the film had reached a value in 
the order of 5 molecular layers, even though an 
intense source of white light was used. Beyond 
this thickness of the film, the sensitivity in- 
creased rapidly. The mercury arc with filters to 
isolate the 3650A band of wave-lengths was then 
substituted for the white light source. 

An accelerating potential of twelve volts was 
found sufficient to saturate the photo-current in 
all cases and was employed in obtaining the data 
for Fig. 2. 

As the thickness of the film increased, the 
photoelectric sensitivity increased rapidly until 
the deposit of sodium reached a value of approxi- 
mately 10 molecular layers. From 10 molecular 
layers to 50 layers, the increase in sensitivity was 
not as rapid as in the earlier deposits. The 
sensitivity then increased more rapidly for thick- 
nesses of film greater than 50 molecular layers, 


until a deposit of approximately 80 layers was 
Beyond 


reached. this, the photo-current re- 
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Fic. 2. Variation of photo-current with the thickness 
sodium film deposited on aluminum. 


mained practically constant for increasing thick- 
ness of film. 

Spectral distribution curves were obtained 
for several films of different thicknesses—from 
fifteen to one hundred molecular layers. Although 
the sensitivity of the thin films was not great 
enough to investigate the curves as far toward 
the long wave-lengths as in the case of the 
thicker films, the results gave no indication of a 
change in threshold (or work function) in going 
from one thickness of film to another. 

A spectral distribution curve for a one hundred 
molecular layer film is shown in Fig. 3. The ordi- 
nates give the photo-current per unit light 
intensity in arbitrary units. The slits of the 
monochromator had to be opened to a con- 
siderable width in order to obtain enough light 
intensity to give an appreciable photoelectric 
effect. This means that a rather wide wave-length 
band (approximately 200A) was used. Fig. 4 
shows how these observations fit the theoretical 
curve of Fowler.’ The points obtained for the 
short wave-lengths fit the theoretical curve very 
well, but those obtained for the long wave- 
lengths show a wide divergence from the theory. 
The horizontal shift of the experimental curves, 
as plotted according to the method outlined by 
Fowler, to make it coincide with the theoretical 
curve, indicated a value of 5150A for the true 
threshold of the sodium film. It was assumed 
that the divergence of the points taken for the 
long wave-lengths from the theoretical curve 


7 Hughes and DuBridge, Photoelectric Phenomena (1932), 
p. 248. 
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was due to the wide band of wave-lengths from 
the monochromator, with perhaps some accom- 
panying scattered light of short wave-length. 
The dotted curve in Fig. 3 was obtained by 
choosing values of the photo-current which would 
cause the longer wave-length points to fit the 
theoretical analysis curve. The contact potential 
between the sodium film and the aluminum 
collecting surface may be obtained from the 
value of the threshold of the sodium film and 
data taken from the current-voltage curves, 
Fig. 5. These curves cut the applied potential 
axis at approximately 0.7 volt. The energy in 
the exciting light (3650A) is 3.38 volts. The 
energy of the fastest moving photoelectron is 
the difference between 3.38 volts and the work 
function of the sodium film. The threshold, 
5150A, corresponds to a work function of 2.40 
volts. The energy of the fastest photoelectron is 
then (3.38—2.40) 0.98 volt. This value when 
added to the 0.7 volt gives the contact difference 
of potential between the two surfaces, or (0.7 
+0.98) 1.68 volts. The threshold of the alu- 
minum was then determined by adding the 
contact potential to the threshold of the sodium 
film, or (1.68+2.40) 4.08 volts. This gives, 
approximately, 3020A as the threshold of alu- 
minum. This value is considerably lower than 
any reported by previous experimenters.* 

Fig. 6 shows the theoretical analysis of the 
data obtained from the current-voltage curves. 
Because the results from the films of different 
thicknesses agree near the “‘tail’’ of the current- 
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Fic. 3. Spectral distribution curve for a 100 molecular 
layer sodium film. The dotted curve represents the curve 
corrected for a wide wave-length band. 





* Reference 7, p. 74. 
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voltage curves, only the results for the thick 
film (100 molecular layers) are shown. When 
the values of the saturation current for the 
different films are reduced to a common level, 
the “‘tails” of the curves coincide for every 
thickness of film. The current-voltage curves 
for two films, as shown in Fig. 5, illustrate this 
point. 

The data for the full line curve in Fig. 6 were 
taken from a paper by Roehr.’ The experimental 
results were plotted in the form log (//x) against 
x as outlined in DuBridge’s theory.'® The hori- 
zontal shift gave a value of 0.96 volt as the 
maximum energy of the photoelectrons at abso- 
lute zero. This compares quite closely with the 
0.98 volt as determined from the analysis accord- 
ing to Fowler’s theory. 

Figs. 7, 8, 9, and 10 indicate the way the 
photoelectric current varies immediately after 
the heat is removed from the sodium reservoir. 
The zero point on the time axis represents the 
instant the photo-current began to decrease after 
the removal of the heater. Figs. 7 and 8 indicate 
that the decrease is about the same for films 
which are less than 50 molecular layers. Those 
greater than 50 layers do not show as rapid a 
decrease in photo-current after a deposit is 
made. It is interesting to note that Andrade and 
Martindale" have presented evidence that the 
films they were studying began to form a definite 
crystal lattice after the thickness had reached 50 
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Fic. 4. Theoretical analysis of the spectral distribution 
curve according to Fowler's theory. 


®W. W. Roehr, Phys. Rev. 44, 866 (1933). 

LL. A. DuBridge, Phys. Rev. 43, 727 (1933). 

" Andrade and Martindale, Phil. Trans. Roy. Soc. 235, 
69 (1935). 
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according to DuBridge’s theory. 


molecular layers. The earlier deposits were aggre- 
gates without a crystalline structure. 

The decrease in photo-current for the early 
deposits is most likely due to the surface motion 
of the deposited sodium atoms. This results in 
the formation of aggregates, which, in turn, lead 
to a reduction in the effective area of the sodium 
which is exposed to the light. 

Bosworth” has shown, in his studies of sodium 
on tungsten, that the decay in the photoelectric 
current is due to a diffusion of the sodium atoms 
into the tungsten. After 10!7 atoms/cm? have 
been deposited the tungsten becomes saturated 
and the photo-current no longer decreases with 
time. Because of the final constant value of the 
photo-current after standing a short time it is 
concluded that no such phenomenon occurs in 
our experiment. 


2 R. C. L. Bosworth, Proc. Roy. Soc. A150, 58 (1935). 
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7. Decrease in photo-current with time after forming a 
deposit of 14 molecular layers. 
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Fic. 8. Decrease in photo-current with time after a deposit 
of 24 molecular layers. 


CONCLUSION 


The essential differences between the photo- 
electric properties of Cs, Rb, or K films on 
silver and the films of sodium on aluminum are: 


1. In the former case, there is a strong photoelectric 
effect when the films are less than one molecular laver 
thick, while in the latter case, the effect cannot be detected 
even with a strong light until the thickness is in the order 
of 5 molecular layers. 

2. The other films show a maximum sensitivity when the 
film thickness is in the order of a few molecular layers 
whereas the sodium films on aluminum do not reach their 
maximum sensitivity until the thickness is in the order of 
80 molecular layers. 

3. After passing through a maximum sensitivity, the 
other films show a decrease in photo-current in going to 
greater film thicknesses. The sodium-on-aluminum films 
show no decrease in photo-current with increase in film 
thickness after the thickness for maximum sensitivity had 


been passed. 
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Fic. 9. Decrease in photo-current with time after a deposit 
of 50 molecular layers. 


The outstanding difference from the theoretical 
standpoint is this—in one case the work function 
of the underlying surface is greater than the 
ionization potential of the alkali metal, while in 
the other case the work function is less. Ives and 
Briggs'® have shown that the electrons originate 
in the underlying metal for the early stages of 
film development in the case of sodium on silver. 
As the thickness of the film increases, more and 
more electrons originate in the alkali metal. 

Because of the comparatively small difference 
between the ionization potential of sodium and 
the work function of silver, there is a certain 
probability that sodium ions will exist on the 
surface of the silver during the early stages of 
deposition. These ions allow the electrons in the 
silver to be released by visible light by reducing 
the work function of the silver. In the case of 
sodium on aluminum, the difference between the 
ionization potential and work function is much 
greater and the probability of adsorbed ions 
being deposited is much less than in the above 
case. The photoelectrons must then originate 
solely in the sodium film. As a result, the sodium 
must be allowed to build up to a relatively thick 


'3 Ives and Briggs, Phys. Rev. 40, 802 (1932). 
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Fic. 10. Decrease in photo-current with time after a de- 
posit of 100 molecular layers. 


deposit in order to present an appreciable volume 
for the absorption of light before a detectable 
photoelectric effect could be expected. 

The evidence from the spectral distribution 
curves and also the current-voltage curves indi- 
cates that the change of sensitivity with the 
thickness of the film is mot due to a change of 
the work function of the surface. This would lead 
one to conclude that the change in sensitivity is 
due to the change in the amount of light absorbed 
in the sodium film as its thickness is changed. 
One need not conclude, however, that the photo- 
electrons in a thick film come from a depth of 80 
molecular layers. The monomolecular layer was 
defined on the basis of a perfectly smooth 
underlying surface. Because of the method used 
in evaporating the aluminum, the underlying 
surface was undoubtedly far from being perfectly 
smooth. The actual surface may be much larger 
than the apparent surface. This point is dis- 
cussed by De Boer." 

The senior author (J. J. B.) wishes to express 
his indebtedness to the Physics Department of 
the University of California for the loan of the 
monochromator used in this research. 


4 J. H, De Boer, Electron Emission and Adsorption 
Phenomena (1935), p. 133. 
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The 1s2s 'S state of He has been treated with a variant of Hylleraas’ variational method, and 
a wave function has been determined which gives an error of only 0.005 ev in the energy of the 


system. The properties of desirable functions for use in this method are discussed. Wave func- 
tions of the Hartree and Fock types, and a more flexible function not involving 712, are also 
determined by simple variational processes. The properties of the resulting functions are com- 
pared, and are discussed in the light of the polarization method of Hylleraas and Bethe, properly 
applicable to more highly excited states. The computations by the polarization method have 
been extended, and the results are compared with those given by the variational procedure. 


INTRODUCTION 


NE of the most interesting of the singly 

excited states of He, as well as one of the 
most difficult to treat, is the 2s state of para- 
helium. The difficulty arises from the fact that 
in the treatment of the moderately excited S 
states of He, methods of approximation which 
are satisfactory elsewhere tend to break down. 
This forces one to rely on variational methods, 
which are in general not so easily applied to 
excited states as to ground states. The best 
previously found wave function for this state is 
that of Hylleraas and Undheim.' The agreement 
of their energy with the experimental value is 
markedly poorer than that given by the Hyl- 
leraas method for the normal or the 2s ortho 
state of the atom.” In the study of this state of 
He, as a preliminary to the treatment of the 
ground state of lithium, we have therefore 
attempted to extend and improve the method by 
increasing its flexibility with respect to the radial 
coordinates, introducing two modifications. The 
first and more important consists of the use of 
two sets of terms with different values of the 
‘outer”’ electron. This 


‘ 


exponent applying to the 
corresponds to the work of Wilson* and others 
who have constructed single orbitals for excited 
electrons by superposing terms with different 
exponentials, finding that only a very few terms 
are required. The second consists in forming the 
wave function as a sum of terms which differ 
somewhat in form from those used by Hylleraas 


E. A. Hylleraas and B. Undheim, Zeits. f. Physik 65, 


1 
759 (1930). 
2 E. A. Hylleraas, Zeits. f. Physik 54, 347 (1929). 
3E. B. Wilson, J. Chem. Phys. 1, 210 (1933). 


and Undheim, which be constructed as 


sums and differences of our terms. Use of these 


may 


terms somewhat mars the elegance of the mathe- 
matical treatment, but contributes to the clarity 
of the physical picture. It also tends to reduce 
the number of terms needed to attain a given 
accuracy in the treatment of the system. 
Having thus produced a considerable increase 
in the accuracy with which this state can be 
treated we have chosen it as one suitable for a 
detailed consideration of the way in which 
atomic wave functions vary with the distance 
between the electrons, and for a comparison of 
the various simpler 
methods of attacking the problem. For this 
purpose it is particularly satisfactory in that the 


successes obtained by 


relation of the electrons is similar to that ob- 
taining in interesting atoms such as Li, Na, 
etc., while this system is much easier to treat 
accurately. 


FoRM OF WAVE FUNCTIONS 


The wave functions employed were of the 
form Y=DomnpCmnpLmnp], where [m, n, p] 
= (x? /4m) fe-"*"! (xr,)™e—9*"2 (Kre)" +e7-"*"2(Kre)” 
e~**"1(xr,)"{(xkrio)”. « is here a scale constant, 6 
and y are parameters, and the distances which 
define a configuration are represented in the 
usual way. Other terms were also introduced, 
which we designate by [m*np], in which the * 
is not an adjunct of the m, but indicates that 
in these terms v has been replaced by another 
parameter, v*. While the best value of « could 
easily be determined during the final variation, 
we desired to base all our computations upon a 
single set of tables computed for fixed values of 
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the parameters 6, v, v*, and after experimenting 
with some simple functions (of the type later 
called Fe), selected the following as suitable: 
6=2, v=0.55, v*=1. Thus, 6 is characteristic of 
a one-quantum orbital about the bare nucleus, 
v corresponds to a two-quantum orbital governed 
by a highly shielded nucleus, and the star terms 
are useful in dealing with the region in the atom 
where the outer electron makes a close approach 
to the nucleus, penetrating a region where the 
shielding is less effective. 

One advantage of the form chosen for terms in our wave 
function appears in the ease with which orbital functions 
can be represented. Thus, such a function as that employed 
by Zener‘ is simply [100]+c[000]. The terms used by 
Hylleraas and Undheim are linear combinations of ours, 
[100 }]+[010], [200}+2[110]+ [020], [200]— [020], etc.; 
it would require three of these to give the Zener function. 
lost in going to more 


This advantage is only partially 


general functions, for the helpful terms are those which are 
not too different from the ones occurring in the product 
functions. Our method enables these terms to be intro- 
duced by themselves, whereas in Hylleraas and Undheim’s 
functions the helpful terms are linked with harmful ones, 
which can be removed only by using additional terms in 
order to cancel them out. For example, in constructing 
our best function, we can do as well by the inclusion of 
[200] and [110] alone as with the last three of the above 
combination terms; the term [020] is harmful and must be 


eliminated. 


The form of trial function adopted enabled us 
to determine analytical approximations to 
several simplified functions as well as to the 
correct function, permitting various interesting 
comparisons. The four functions which our 
approximations approach as limits are: 

The Hartree function (//#); the best function 
expressible in the form K(r,) X L(re), where K(r) 
and L(r) are single-electron functions of the 
general character of 1s and 2s orbitals. 

The Fock function (F#); the best function 
which can be expressed as K(r;) X L(re)+K(re) 
xL(r:), where K(r) and L(r) are similar to but 
not identical with those occurring in the //®. 

The Best Radial (R®); the 
function of any form containing only 7; and fo, 
which must enter symmetrically. 

The Correct function (C®); a function of 1, 
re, and rjo, the first two entering symmetrically. 

(The above abbreviations will be applied also 


function best 


*C. Zener, Phys. Rev. 36, 51 (1930). 


TIONS 
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to our analytical approximations where no con- 
fusion will result.) 

It must be clearly understood that these four 
limiting functions have a definite existence inde- 
pendent of the method used to approximate 
them. Thus R® is the exact solution of a differ- 
ential equation obtained from the Schrédinger 
equation upon replacing 1/r;2 at each point in 
configuration-space by its mean value over all 
points having the same values of r; and fo, i.e., 
by 1/g:2 where ge is the greater of r; and fo. 
There results an equation involving only r; and 
ro, but it cannot be separated and solved by a 
product of single-electron functions. The equa- 
tions solved by F and J//# are sufficiently well 
known. 

In the ensuing discussion we shall distinguish 
between two forms of atomic wave function 
constructed from the Hartree single orbitals. 
The first, which we shall designate //)?, is that 
already described, which neglects the require- 
ment of symmetry in electron coordinates. The 
second, //,4, meets this requirement by the 
addition of exchange terms: it has thus the same 
form as F4, but the coefficients are those deter- 
mined by minimizing the energy of // 9. To find 
the energy of //, requires an additional com- 
putation, using these coefficients explicitly. 

In using the variational method to find ap- 
proximations for the various types of limiting 
function, we have simply to observe certain 
restrictions on the character of the terms used 
or the way in which they enter. The product 
functions 7 and F# must have the form 


Dd md n€mn' Lm, 2, OJ + DY md nCm*en'[m*, n, 0) 


corresponding to K(r)=e-*’Soc,’r", L(r)=e-"" 
YLenr™ +e" Yc, *r", (taking «=1). In deter- 
mining //#, of course, the exchange terms have 
to be omitted from the definitions of the func- 
tions and from the matrix formulas. Correspond- 
ing to the conventional way of determining these 
functions, we first assume that K(r)=e-*’, the 
He* 1s function; we thus determine preliminary 
values of the c and c* by variation among all 
terms having »=0, p=0. The next stage of 
approximation improves K(r), determining the 
c’ by a variation among individual terms each 
in the form >»fcn{lm, n, 0]+c,*[m*, n, 07}. 
The c and c* can then be readjusted if necessary, 
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holding the c’ constant, and so on until self- 
consistent results are reached. R® is obtained by 
an unrestricted variation among all terms having 
p=0; the removal of this specialization gives C®. 

In none of these variations (except that deter- 
mining the K(r) factor in the product functions 
H# and Fé) were we seeking an absolute mini- 
mum value of the energy parameter. Hylleraas 
and Undheim have shown, however, that the 
second root of the secular equation can never lie 
lower than the second characteristic energy value 
of the problem, and approaches that value as a 
limit as the terms in the series are improved, even 
though the lowest root may be at the same time 
a very poor approximation to the lowest energy 
value. The only precaution in using the method 
is, therefore, to be sure that one has the proper 
root. 


INTRODUCTION OF 7j2 


There has been considerable experimentation 
to determine the most practical way of intro- 
ducing fiz into the wave function. It will not do, 
(as Hylleraas and Undheim point out) merely to 
multiply a function of the R® type by (1+ cri). 
Inspection of Table I will show that we find no 
tendency for the coefficients of the terms 
[m, n, 1] to be proportional to those of the cor- 
responding terms [m,n,0], a result likewise 
noticed by Hylleraas and Undheim. A possible 
guide to the selection of 712 terms can be derived 
by applying a method due to Hassé and Lennard- 
TABLE I. Computed wave functions and energies. 


FUNCTION Hot H,® Fe R& Ce 








Energy, in 2Rh —2.13915 —2.14309 —2.14407 —2.14580 

Energy error (ev) 0.185 0.078 0.052 0.005 
1.00 1.00 1.05 1.05 
4.223 16.756 16.485 16.036 


5 perelectron 14.447 1 


COEFFICIENTS OF TERMS IN NORMALIZED FUNCTIONS 


(000) 25.21939 8.41368 4.88925 

(100) —6.32769 -—2.75611 —0.02948 0.35368 
(010) —0.00725 —.85415 — 69387 
(110) .00237 .23488 . 14042 
(200) 0.28620 .05770 .09866 -12251 
(020) .03635 03271 

(210) — .00005 

(120) —.01191 

(220) .00025 

(0*00) —20.80546 —5.66385 —7.17757 —2.30134 
(1*00) —7.90329 —2.14445  —2.26128 0.47626 
(0* 10) .00488 1.39627 1.08184 
(1*10) .00185 

(0*20) — .02447 

(1*20) — .00926 

(001) 0.63432 
(101) — .08422 
(011) 07633 
(0*01) — 1.21923 
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Jones.® We regard R® as the solution of an unper- 
turbed equation, and introduce as a perturbation 
the quantity 1/ri2.—1/gi2. Then the Hassé 
approximation to C¢ is 


RO{1+c[1/ri2—1/g12]}, 


c being a constant to be determined by variation. 
The energy corresponding to this function is not 
integrable over small ry», but this is a defect in 
the approximation and not a property of the 
function which we desire to represent. Upon 
multiplication of the perturbation term by 
rio/212 We remove this defect and obtain an 
expression which is but little different from the 
Hassé approximation in regions where one r is 
much greater than the other, these regions being 
those where the approximation itself is valid, 
and simultaneously those in which the factor R® 
takes on its largest values. In a wave function for 
an excited atom, it is a small approximation to 
replace gi2 by r; and re, respectively, in the terms 
in which electron 1 and electron 2 are assigned 
the more diffuse distribution. From each term 
[m,n,0] of Rb we thus derive [m, n, 0] 
+c{[m—1, n,0]—[m—2, n, 1]}. The argument 
not being rigorous, we can in general get the best 
results by introducing the new terms with free 
coefficients to be determined by variation, rather 
than adhering to the relations between the coef- 
ficients which the Hassé approximation pre- 
scribes. Nevertheless we should expect definite 
cooperation between the terms which in the 
theory occur in pairs; each term in 732 may be 
most helpful only when a corresponding term, 
which behaves in nearly the same way for large r, 
is included in the series to be varied. 

Terms suggested by this reasoning cannot be 
used for m<2, for they would not be properly 
integrable. It may therefore be necessary to 
employ, say, the terms [ m,n, 0] and [m—1, n, 1]. 
Moreover, terms such as [m+1*,m,0] and 
[m*,n,1], are always available, and may to a 
certain extent take the place of [m—1, n, 0] and 
[m—2,n,1], since both sets have the common 
feature of representing more compact distribu- 
tions than the term [m, n, 0] from which they 
are derived. We believe that this explains the 
fact that the introduction of star terms is more 





5 J]. E. Lennard-Jones, Proc. Roy. Soc. A129, 598 (1930); 
H. R. Hassé, Proc. Camb. Phil. Soc. 26, 542 (1930). 

















WAVE FUNCT 
helpful in constructing the function with 7). than 
when the terms with r,. are omitted. 

The importance of a given term is by no means 
indicated by the magnitude of its coefficient in 
the series, partly because the terms as we define 
them are very differently normalized, and partly 
because it may happen that two fairly similar 
terms enter with large coefficients of opposite 
sign, the net result being a very small alteration 
of the function which might be ignored without 
serious effects. (In such a case it is not worth 
while to consider the corresponding terms in 112.) 
A principal and indispensable term (as indicated 
by its effect on the computed energy) may thus 
have a relatively small coefficient, and serve as 
the basis for the construction of useful 712. terms, 
while other terms with larger coefficients are 
ignored for this purpose. 


EVALUATION OF MATRIX COMPONENTS 


The required matrix elements can all be found 
in terms of the integrals 
J(f, g, h/a, B)=(1/162?) fd Vrife-2"1re%-F 270". 
2 n+i(a@, B) + Vnsze, moilB, @), 
-2(a@, B) + Vase, me2(B, a) 
= A se(a) XAxni2(8), 


3(B, a) 


Jim, n, —1)= Vas 


J(m, 2, 0) = Vanes. ; 


J(m, n, 1) = Ving2, n+3(a@, B) + Varo, m4 
+30 Vines n41(@, B)+ Vag 


4, m+i(B, a) ], 
J (me, 8, 2) Vass. a4 


2 (a, B)+ Vase. m+4(B, a) 
a fe na2(a, B)+ Vas 1, m+2(B, a) 


=A wis(a)XAns (8B) +Amss(a) XA n+2((). 


Here A,(a)= fo"e-™ X"dd 


and V,,(a, 8)= [ h"e oan pre P"du. 
ei) ei 


The latter quantities are conveniently computed 
by using formulas given by the writers.® 
Consider two terms [m, n, p] and [m’, n’, p’}; 
as either or both may carry stars we must dis- 
tinguish between the possibly different values » 
and v’ which occur, whereas «x and 6 must be the 
same in both terms. Between these terms the 








®H. M. James and A. S. Coolidge, Phys. Rev. 49, 688 
(1936). 
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matrix element of unity is 


2J(m+m',n+n’', p+p'/v+v’, 28) 
+2)(m+n',n+m', p+p'/vt+é, v' +4); 


which may be abbreviated to 2/(000) + 2/’(000), 
the general scheme of abbreviation being such 
that J(a, b, c) stands for J(m+m’'+a, n+n'+), 
pb+p'+c/v+y’, 25) and J’(a,6,c) is similarly 
formed from the “exchange” term. The matrix 
element of the energy, in units of 2R,h, is the 
sum of two terms involving J and J’, respectively, 
with identical coefficients as follows: 
x? { —(v? +8) J(000) + »(2m+ p+ 2)J(—100) 
—m(m+1+p)J(—200)+6(2n+ p+ 2)J(0—10) 
—n(n+1+p)J(0—20)+p[{ —2(p+1)-—m—n} 
xX J(00—2)+»J(10—2) —»J(—12—2) 
+mJ(—22—2)+6/(01 —2) —6/(2—1—2) 
+nJ(2—2—2)]}+2«{J(00—1) —2/(0—10) 
—2J(—100)}. 


It will be noted that these coefficients depend 
only on the first of the two combining terms. A 
symmetrical formula could of course be obtained; 
but the unsymmetrical form has the advantage 
of permitting a check against the alternative cal- 
culation in which the terms are interchanged. 


RESULTS 

In Table I will be found the wave functions 
and energies obtained in our best computations 
with functions of the various types. 

The energies stated are those for the fixed 
nucleus problem, in units of 2R,4=27.074 ev. 
Hylleraas* has shown that, to the approximation 
in which one can neglect the dependence of the 
wave functions on fiz, these numbers also give 
the energies of the corresponding approximate 
solutions of the free nucleus problem in units of 
2Ruyuet=27.070 ev. Bethe’? has shown that for 
the ground state the correction for dependence 
on ry of the function is a fifth of that above. In 
the present computations this further correction 
has been neglected, for the dependence of C® on 
ri2 is much less important; the error in the energy 
arising from neglect of r;. terms is less than a 
tenth of that found for the ground state. The 
errors given in Table I are then differences 


between the corrected energies and the observed 
energy of the system. 


7H. Bethe, Handbuch der Physik, Vol. 24/1. 
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A discussion of these results falls naturally into 
two parts; in the first we are concerned with the 
reliability of the approximations which the vari- 
ation method yields, taking the four limiting 
functions and their corresponding energies as 
standards of perfection, while in the second it is 
these limiting functions themselves which inter- 
est us, and which we regard as sufficiently well 
represented by our approximations. 


ACCURACY OF RESULTS OBTAINED BY THE 
VARIATIONAL PROCEDURE 

Hylleraas and Undheim estimate that the 
limiting energy obtainable by a function of the 
radial type is — 2.14400. We have barely sur- 
passed this, reaching — 2.14407. There appears 
to be no calculation of the function itself avail- 
able for comparison. This is also true of the 
correct function. 

As measured against experiment, the error in 
our approximation to the energy of C® is only 
about one-sixth of that associated with the six- 
term function of Hylleraas and Undheim, the 
best previously available. The essential difference 
between the two functions is the addition of the 
v* terms, to be discussed later. Other terms were 
introduced, but discarded as being of no value. 
The explorations made were, however, not suf- 
ficiently extensive to permit an accurate estimate 
of the convergence limit, in accordance with the 
procedure outlined elsewhere. * 

In the case of the Hartree function, it is 
possible to make a direct comparison between the 
results of the variation method and those ob- 
tained by accurate numerical integration by 
Wilson and Lindsay® at the first stage of the 
calculation, in which the 1s orbital is taken as 
e~*". The energy parameter ¢ for the 2s orbital 
we find to be «=0.3054, and Wilson informs us 
that they obtained precisely the same value, this 
proving the effectiveness of our variational pro- 
cedure. In view of the existence of their complete 
solution, we have not carried out any further 
stages of approximation. In any case, the dif- 
ferences between our function and their final 
one are so small that we shall use our function 
in the following discussions. 


*H. M. James, A. S. Coolidge and R. D. Present, J. 


Chem. Phys. 4, 187 (1936). 
®W.S. Wilson and R. B. Lindsay, Phys. Rev. 47, 681 


(1935). 
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The value given by Wilson" for the energy of 
this state is, however, very different from ours, 
his computed ionization energy being less than 
that observed by 0.81 ev, while ours is in error 
by only 0.184 ev. Wilson calculated the energy 
of the atom by introducing into the formulas of 
Slater" values which he obtained for the integrals 
using his Hartree orbitals. In doing so he neg- 
lected the nonorthogonality of his orbitals; his 
value is accordingly not the energy associated 
constructed the 


with a wave function from 
Hartree orbitals, as ours is, and his computed 
ionization energy is not comparable with other 
results of the Hartree method. 

Our approximate Fock function was obtained 
as a second approximation, varying each factor 
once. A further variation gave no appreciable 
change in the 2s orbital; this function thus 
represents the best approximation to the Fock 
function obtainable with the terms included. It 
may be compared with that computed by Hyl- 
leraas,” likewise by the variation method, but 
with a different series, and without varying the 
1s orbital. Taking »=0.5, he finds with the three 
terms [m00], m=0, 1, 2, the energy — 2.14275, 
which sinks to — 2.14294 upon inclusion of two 
more terms of the same series. Judging by the 
smallness of the improvement given by the fifth 
term, he accepts this value as substantially the 
limit of convergence. With the same three terms 
as Hylleraas (but with »=0.55) we find the 
energy —2.14285; putting [0*00] in place of 
[200] gives us —2.14293, illustrating the use- 
fulness of the star terms. With [200] and [1*00] 
in addition, we reached the value — 2.14307, 
appreciably lower than Hylleraas’ apparent 
convergence limit. The value given in Table I, 
— 2.14309, shows the minute improvement re- 
sulting from variation of the 1s orbital, and 
justifies Hylleraas in neglecting this step. Com- 
parison of numerical values of his 2s orbital with 
ours shows typical discrepancies of one or two 
percent; in view of the lower energy, we believe 
our result is superior. Our complete atomic 
function should be more reliable than that of 
Hylleraas also because we have improved the 1s 
factor. 


10 W.S. Wilson, Phys. Rev. 48, 536 (1935). 
17. C. Slater, Phys. Rev. 34, 1293 (1929). 
2 E,. A. Hylleraas, Zeits. f. Physik 83, 739 (1933). 
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COMPARATIVE ACCURACY OF THE METHODS 
OF APPROXIMATION 


The functions //,6 to C@ which we have 
obtained may be considered as successive steps 
in the determination of an accurate function, in 
which approximations are successively removed. 
The approximations in //, consist of neglect of 
exchange in determining the function and failure 
to give a detailed treatment of the interelectronic 
repulsion. In Fé the first of these is removed, 
while R® and C# represent two steps in taking 
better account of the electronic interaction. For 
comparison with our results on this state there 
appear in Table II certain corresponding values 
for the ground states of He and Be. 

The exchange effects in 152s \S He and (1s)?(2s)? 
Be are comparable with regard to the quantum 
numbers of the states involved; in Be, however, 
there are twice as many exchange terms. That 
the exchange effects in these two states differ by 
roughly a factor of 2 in absolute magnitude is 
thus not surprising. The percent of error in //,® 
due to the neglect of exchange is, however, very 
different from that in Be, where there is not only 
an error due to neglect of the polarizing effect 
of the outer electron on the core, which has its 
analog in the excited He, but also a very large 
error due to failure to treat the interactions of 
the 2s electrons in a detailed manner. In such a 
case as this the Fock method does not offer a sig- 
nificant improvement over the Hartree method. 

Abandonment of the orbital idea, the difference 
between Fé and R®, permits a partial adjust- 
ment of the relative positions of the electrons. 
In both of the states of He a third of the error 
of the Fock function is thus removed. This may 
be taken as a rough measure of the general 
effectiveness of this modification. Thus, for 
instance, a function similar to R® applied to the 
treatment of Be might give an error in the first 
two ionization energies as small as 1.0 ev. 
Further reduction of this error, associated prin- 
cipally with the treatment of the 2s electrons, 


TABLE II. Errors in computed energy necessary for removal 
of two electrons. 


1 S He is2s'S He 1 2 Be 
H,® 10 0.185 1.95 
Fo 1.14 0.078 1.61 
Re 0.70" 0.052 


would require the very difficult step of intro- 
ducing explicitly the distance between the 2s 


electrons. 
RADIAL CHARGE DENSITIES 


It is something of a problem to find a satis- 
factory method for portraying the wave functions 
of the different types, with the special object of 
facilitating intercomparisons. Fig. 1 shows the 
total radial charge densities; in order to avoid 
confusion due to overlapping, we have plotted 
a single curve for our best function, and indicated 
by an enlarged difference plot how the other 
values deviate from this. C@ gives a charge 
density of the form p,=?r’>-,, ,C,,e~%«'r’. The 
coefficients are shown in Table III]. Though 
useful for certain purposes, the total density is 
not adequate for a careful study of the functions, 
since it suppresses the most interesting charac- 
teristics of the better functions, which take into 
account in detail the interaction of the electrons. 
Thus, although the total density given by R® 
is much better than any of the others without 
’\o, it would be a mistake to conclude that R® 
is itself a correspondingly close approximation 
to the correct function, as consideration of the 
errors in the energy will show. On the other hand, 
in cases where only the charge density is of 
interest it may be a helpful simplification to use 
functions of this type. Improvement of the 
dependence of the functions on interelectronic 
distance will be seen to result in an increase in 
the charge density for the smaller values of r, 
when the electronic repulsions are most im- 
portant. This is probably a general behavior. 


a ae n — H 
’ ‘ x He 
© 091 4 . . H 


Fic. 1. Total radial charge densities. Upper figure shows 
values for Cé. Lower figure shows values for other func- 


tions minus those for C&@. Normalized to 2 on integration 


over r. 
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raBLe III. Constants in formula for total charge density of C® pp=rE¢,7e7°o"Cop?’. 


Cc 


g 5g og 1 ao 

1 1.155 -0.02237 0.07 185 
2 1.6275 .01346 01260 
3 2.1 .03986 23617 
4 2.6775 .18754 —.12678 
5 3.15 —.19766 70711 
6 4.2 — .06098 37.37703 
7 5.355 — 12.83385 —8.66312 
8 5.8275 14.76190 6.808 16 
9 6.3 — 1.88790 1.48937 

On the other hand it cannot be expected in 


general that, as is observed here, the true radial 
density will be between those given by H® and 
Fé, for as we shall see, the relation between 
these functions varies with the system treated. 


COMPARISON OF HARTREE AND FOCK 
FUNCTIONS 


The product functions H@ and F# can be 
compared by plotting or tabulating their separate 
1s and 2s factors. Table IV gives the results for 
Fe, the quantities tabulated being the “radial 
amplitudes” c\rK(r) and cerL(r), with c; and ce 
so chosen that each factor separately will be 
normalized with respect to integration over r 
alone. Since these two factors are not orthogonal, 
these normalization constants are not the proper 
ones for constructing a normalized atomic func- 
tion; for the same reason, caution should be 
observed in comparing these orbitals with those 
given by other authors who have used the Fock 
method in problems where the symmetry was 
such as to permit solution in terms of orthogonal 
orbitals. In such cases comparison should be 
based on the complete atomic functions. Fig. 2 
shows our 2s factors of //@ and F4; the 1s factors 
would on this scale be imperceptibly different 
from the Het orbital, and are not shown. The 
figure shows that the 2s factor of F® is more 
diffuse than that of //#; this difference of course 


TABLE IV. Radial functions in Fe, 

r ewK(r) corL(r) r cwK(r) carL(r) r crL(r) 
0.2 0.756 0.1033|2.4 0.114 -0.2487] 64 -—0.3787 
04 1.014 0.1513} 2.6 0.083  -0.2890| 7.2 0.3205 
0.6 1.021 0.1601/2.8 0.061 -—0.3261/ 8.0  —0.2615 
08 0.913 0.1422/ 3.0 0.043 -0.3582| 9.0 —0.1949 
1.0 0.766 0.1066 | 3.4 0.022 -0.4071] 10.0 —0.1406 
1.2 0.617 0.0608 | 3.8 0.011 —0.4406/ 11.0 —0.0987 
1.4 0.484 0.0086} 4.2 0.005 —0.4573| 12.0 —0.0671 
1.6 0.372 -0.0459 4.6 0.003  -—0.4602/ 13.0  —0.0456 
1.8 0.281 —0.1006 | 5.0 —0.4533 | 14.0 —0.0302 
2.0 0.210 —0.1528 5.4 —0.4388 | 15.0 —0.0197 
2.2 0.155 —0.2026 | 5.8 —0.4181 } 


( v1 o2 o3 o4 
0.10055 0.10313 0.007 38 0.000 13 
—.19053 —.14547 — .00476 

. 16818 .04403 

—.16386 —.10751 .05943 00199 

.19377 —.01275 — .03667 

2.90319 1.26657 02744 02464 
— 2.91047 —.57610 01630 

1.54834 .72437 

1.68280 


carries over into the total charge density, since 
that due to the 1s factor and to the exchange is 
nearly the same in both cases. Fig. 1 shows that 
the total charge density due to F® is somewhat 
excessive at large r, while both 77)? and //,% err 
more seriously in the opposite direction. These 
results contrast with those of Hartree and 
Hartree,'* who find for Be that the Fock density 
is more compact than that given by the self- 
consistent field theory, and surmise that this is 
a general relation. The origin of the difference 
is not far to seek. In the case of excited He, the 
exchange term occurs with a positive sign, just 
as in the Heitler-London treatment of two 
attracting H atoms; and for the same reason: 
there exists a degeneracy of configurations of the 
same total spin. In Be, and in general in atoms 
where the inner shells are filled, the situation 
resembles that of two repelling He atoms, in 
which there is no degeneracy of configurations 
and the exchange terms have negative coef- 
ficients. It is not surprising that the effects of 
considering exchange should be of opposite 
character in these contrasting cases. Similar 
remarks apply to the mean values of f° cor- 
responding to Fé and H,®. We find the former 
some 16 percent greater (see Table I), while in 
Be it is 11 percent smaller. (The Hartrees give 
8.42 and 9.54 as the respective contributions 
from each electron in the 2s orbit. A rough 
numerical integration of their 1s function gives 
0.23 as the corresponding 7*, making the mean 
values per electron 4.32 for F& and 4.88 for 
H,®.) This result, however, is to be expected 
only in excited atoms, and therefore is hardly 
an objection to the general usefulness of the 
statement made by the Hartrees, that wave 
functions calculated by the method of the self- 


18D. R. Hartree and W. Hartree, Proc. Roy. Soc. A150, 
9 (1935). 














WAVE FUNCTIONS FOR HELIUM 683 











Fic. 2, 2s orbitals of Hartree function (H) and Fock 
function (F). Curves show rL(r), normalized to 1 on 
integration over r. 


consistent field appear always to give too large 
a value for the diamagnetic susceptibility. 


THE POLARIZATION METHOD 


The relation between the Fock function and 
those which are not restricted to the orbital 
form may be qualitatively understood with the 
aid of the polarization theory developed by 
Hylleraas“ and Bethe.’ The method, primarily 
applicable to highly excited atoms, makes use 
of a separation of variables similar to that gener- 
ally used in molecular problems. The outer 
electron is regarded as moving so slowly in com- 
parison with the inner that the motion of the 
latter is essentially that of a particle about two 
fixed centers, and may be represented by a wave 
function w#2(1), a solution of the differential 
equation 


(-—3V"°-Z r;—-Z rot+1, 1 12)Uo(1) = Fous(1 a 


which is regarded as a function of the coordinates 
of electron 1 only, those of electron 2 entering 
as parameters. This function will differ from 
K(1), the Hartree inner orbital, by being shifted 
in a direction away from the instantaneous 
position of the outer electron, so that its re- 
pulsive effect upon the latter is reduced. The 
outer electron is in turn considered as moving in 
the field of the thus polarized inner distribution 
plus that of the nucleus; its motion will be de- 
scribed by (2), a solution of a non-Coulomb 
central field problem. Clearly v(2) will be more 


4“ E. A. Hylleraas, Zeits. f. Physik 66, 453 (1930). 


compact than the corresponding Hartree orbital, 
in determining which the polarization of the 
inner orbital was neglected. As an approximate 
solution of the wave equation one then takes the 
definite function ~=wue(1)v(2); there are no 
further steps of successive approximation, so 
long as polarization is dealt with in this manner. 

In the case of a moderately excited atom the 
reasoning loses sharpness. In the first place, 
there is a considerable probability of the electrons 
being at comparable distances, so that neither 
can be regarded as almost stationary. In such 
configurations, we must (in terms of the classical 
analogy) think of the orbit of each electron as 
somewhat polarized by the instantaneous prox- 
imity of the other, but to only a fraction of the 
extent which would be caused by a stationary 
electron, since the actual electron does not 
occupy any given position long enough for the 
full effect to develop. The polarization correc- 
tion as given by the Hylleraas-Bethe theory will 
therefore be excessive in the region of penetra- 
tion. In the second place, the exchange terms 
required by symmetry are not negligible, and 
their interpretation is not clear. It is assumed 
that Y= we(1)v’(2)+,(2)0’(1), in which w(1) is 
the function obtained above, while v’(2) is the 
solution of the variational problem in which the 
energy of the complete function is minimized. 
To a first approximation the corrections due to 
exchange will be independent of the polarization 
corrections, and one may expect v(2) and v’(2) 
to be related in the same way as the outer orbitals 
given by the Hartree and Fock methods, 
respectively. 

We shall now show that even in the compar- 
atively little excited state which we are now con- 
sidering, the qualitative, and to some extent the 
quantitative, properties of our functions are 
those indicated by this theory. 

An exact calculation of u2(1) would be very 
difficult. Bethe’? and Ludwig" have applied a per- 
turbation method, starting with the simple Het 
function w4°(1), and calculating the first-order 
perturbation produced as the outer electron 
moves in from infinity. They use an expansion 
in surface harmonics, and obtain only the terms 
in P; and P», the first being much the more im- 
portant. We have obtained a general formula for 








6 (4, Ludwig, Helv. Phys. Acta 7, 273 (1934). 
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these terms, and have included in our computa- 
tations terms in Py and P3. (Details are reserved 
for a separate section.) The new terms are small 
except when the second electron is only slightly 
beyond the first, or even inside it,—that is, for 
configurations for which the basic ideas of the 
method break down anyway. 

The function v(2) has received practically no 
discussion by previous authors, who have been 
concerned primarily with determining energy 
values, for which purpose they have used Fock 
or even hydrogenic outer orbitals, as being suf- 
ficiently good in the case of highly excited states. 
We have not attempted to determine (2) 
directly, but have investigated an ‘“‘empirical”’ 
v(2) of form so chosen that, when combined with 
our approximate values of u2(1), it reproduced 
C® as well as possible. In the outer regions, this 
function is well defined by the total radial 
density; since the exchange terms are negligible 
here and u2(1) is normalized for integration over 
dr, the density gives immediately the square of 
v(2). In the inner region, the choice of 2(2) 
becomes increasingly arbitrary; no values can 
be found which give consistent agreement, and 
those selected will serve as well as any to show 
the magnitude and nature of the discrepancies. 

In Table V are indicated values of the assumed 
v(2), the corresponding atomic function P®, and 
Cé, for several configurations. We also include 
the angle variations (differences in the function 
for configurations with the same r; and fs, but 
minimum and maximum values of 7,2, compared 
with the mean value over all values of 7,2). This 
table shows that, as expected, the function (2), 
in the region where it is well defined, is appreci- 
ably more compact than the Fock outer orbital, 
and that it leads to an atomic function in quite 
satisfactory agreement with the correct function, 
both in absolute magnitude and in angle vari- 
ation. The angle variation furnishes a particu- 
larly good test of the polarization theory, since 
the values calculated are only slightly dependent 
upon the form assumed for v(2). In the inner 
region, on the other hand, the agreement becomes 
much worse, and the angle variation according to 
the polarization computation is much too large, 


likewise as expected. 
Reasoning similar to that above can be applied 
to the determination of the function R&. The 
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inner function #o(1) now describes the motion of 
the first electron in the field of a second electron 
distributed over a sphere with radius re; the 
first-order correction is given by the term in P» 
of the complete polarization theory. The “‘polar- 
ization”’ takes the form of an increasing diffuse- 
ness of u%2(1) as r2 becomes smaller. It must also 
result in some contraction of the outer function 
v(2), compared with Z(2), but the smallness the 
P, term in the polarization correction indicates 
that this should be insignificant, except for very 
small ro. 


COMPARISON OF THE FoOCK FUNCTION WITH 
More ACCURATE FORMS 


Except for the complication introduced by the 
exchange terms, we have now a basis for com- 
parison between R® and F®. In Fig. 3, we have 
plotted 7;r2R® in the form of a series of curves 
for fixed values of 7;, which we may regard as 
cross sections of the surface representing 7\r2R® 
as a function of 7; and re. These cross sections 
are practically indistinguishable on this scale 
from those corresponding to F. For small 7, 
they have essentially the character of a 2s 
orbital, while, as 7; increases, they pass over into 
the 1s orbital form. The exchange effect, a neces- 
sary consequence of the symmetry of the 
function, largely masks that due to polarization. 
It happens, however, that, for r;=0.2, F® gives 
r,K(1) =0.756, 7,L(1) =0.1033, while for 7;=1.0 
we find 7,K(1)=0.766, r,L(1) =0.1066. A com- 
parison of the corresponding cross sections for 
R# and F® therefore practically eliminates the 
exchange effect, since the amounts of 1s and 2s 
orbitals present in F are almost the same in 
both cases, and the exchange effect is supposed 
independent of the polarization correction. These 
cross sections are shown in Fig. 4. For large fe 
and small 7;, both F@ and R® reduce practically 
to the single terms K(1)L(2) and wo(1)v(2). 
Further, (1) is, as shown by computation, 
practically indistinguishable from K(1), so that 
in this region we are essentially comparing 2(2) 
with L(2); it is seen that the difference between 
R# and F¢ is practically the same for both values 
of r; and corresponds to a slight contraction in 
the outer function. On the other hand, when 
both r’s are small the two terms in the function 
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Fic. 3. Values of best function of 7; and rz. Curves give r:72R® as functions of rz for the indicated 
values of r;. On the scale shown, R® is normalized to 167? on integration over coordinate space. 








Ar 





Fig. 4. Comparison of R& and F#. Curve F shows values of r2F4, for 7, =0.2 and for r;= 1.0, the 
values being too close together for separate representation. The other curves show, ten times en- 
larged, the values of r2.(R@— Fe) for r; as indicated. Normalization as in Fig. 3. 


are of comparable importance, and both must be 
somewhat affected by polarization. We might rep- 
resent R® as K(1)L(2)g(1, 2)+L(1)K(2)A(1, 2), 
where the factors g(1,2) and h(1, 2)=g(2, 1) 
have, both separately and in combination, the 
character of moving the average distribution of 
electron 2 in or out, respectively, accordingly as 
the known position of electron 1 moves out or in. 
In the Fock function Z(2) and K(2) are dis- 
tributions in the field of the average distribution 
of electron 1. In R&, however, L(2)g(1, 2), as a 
function of 2, is more diffuse than this if the 
fixed electron (considered as a K electron) is at 
an unusually small r(0.2) and less diffuse if it is 
at a greater than average r(1.0). On the other 
hand, K(2)h(1, 2) is in both cases more diffuse 
than the Fock orbital K(2), since 0.2 and 1.0 
are both much smaller than the average distance 
of an L electron. The term K(1)Z(2)g(1, 2) is the 
principal part of the function for intermediate 
values of r2, and it may therefore be expected 
that the Fock cross section should lie between 





the other two in this region. For still smaller 
distances the polarization idea breaks down com- 
pletely and the theory ceases to be helpful; it 
remains true, however, that the Fock cross 
section lies between the other two. 

C# differs from Fé, to an even greater extent 
than does R®, in assigning a larger probability 
to the appearance of the two electrons at approx- 
imately the same r. This is, of course, not cor- 
related with a corresponding increase in the 
indicated chance of the two electrons approach- 
ing each other, for the favored configurations 
place the electrons on opposite sides of the 
nucleus. At large r the charge density is distinctly 
more compact than for Fé and R®, as indicated 
by the polarization theory. In our diagram their 
difference in form between R® and C4, really 
more pronounced than that between F¢ and R®, 
is somewhat disguised by the different normal- 
ization of the outer parts of the distributions. 

Table V will illustrate the importance of the 
dependence of C® on rj, particularly for the 
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smaller values of r; and ro. Here the differences 
in C@ with varying r;2. are much greater than the 
differences between R® or F@, and Cé, the mean 
value of C® for given r; and re. R® is in fact a 
fairly good approximation to C& everywhere; 
from this the satisfactory agreement of the charge 
densities for R@ and C# follows. 


DETERMINATION OF FIRST-ORDER POLARIZA- 
TION OF THE CORE 


COMPLETE 


We require a solution of the equation 
(QV 2+E (re) +Z/171+Z/re—1/riz)u2(1) =0, (1) 


where 72 is considered as a fixed parameter, and u2(1) is a 
function of the coordinates of electron 1 which depends 
upon the position of electron 2. We regard (1/r:2—1/r2) as 
a perturbation term, small of the first order. Then, in terms 
of the new variables 


p=Zr,, R=Zn, 


we find the zeroth approximations for solution and energy 
to be 


u%™(1)=(Z/2)ce", c&=4Z/r; E(R)=—}3Z*?9—Z(Z—1)/R. 


To the first order, the energy is E°(R)+«,(R); where 
e:(R) = —Ze-**(1+1/R). The first approximation for the 
solution is assumed to have the form 


@ 
uo(1)=u%(1){1+ = f,(1, 2)P.(cos @12)}. (2) 
n=VU 


To determine the f,(1, 2), we expand 1/rj,. in spherical 
harmonics, substitute (2) in (1), and equate terms of cor- 
responding order of smallness and involving the same 
harmonics. Let There result the 


equations 


wu (1)f,(1, 2)=¢,(p). 


go — got2(¢o' + ¢0)/p 
pf e **(1+1/R) for p<R, (3) 
=Cce ° 
. L[e?*(1+1/R)+1/p—1/R] for p>R. 


On’ — On t2(en'’ +en)/p—n(n+1)/p? 
ae" R"*, p<R,\ " 
LR*/p"*, p>R,J 


Ludwig gives the solutions of these equations for n=1, 2. 
(Through an error, he introduces a factor 2 on the right 
side, but his solutions are correct except that in his Eq. 
(15a) (1+ R?) should read (1+ R)?.) The general solution for 
n>0O may be obtained in the same way. 

The case »=0 presents some interesting peculiarities. 
For the region p<R, a particular solution, everywhere 


=ce >O. (4) 


finite, is found to be 
Hin = he *®(1 +1/R) ce? Ei(2p) 

—log (2p)—p—(e”—1)/2p]}. (5) 
To this may be added a constant, a(R), times that solution 
of the homogeneous equation which remains finite for 
p—0. Similarly, for p>R a particular solution, vanishing 
properly at infinity, is 


Yout = }ce Pile 22(1+1/R)+1-1 R} 
X[1/2p—log (2p)—p]+p}. (6) 
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FOR HELIUM 
To this may be added a constant, 8(2), times that solution 
of the homogeneous equation which vanishes properly for 
p—«. One would expect the constants a and £8 to be 
determined by the conditions that the two solutions and 
their derivatives must have the same values at p=R. 
However, it turns out that the two required solutions of 
the homogeneous equation are identical, namely, «(1 


itself. Continuity of function requires 


B—a=f°=(1/Z){(e*®Ei(2R)—1)(14+1/R 
+log (2R)(1—1/R)}, (7) 


and it is found that this relation also leads to continuity of 
derivative. 

There are therefore an infinite number of solutions of the 
complete equation which satisfy the boundary conditions 
The one desired is that which is orthogonal to the unper- 
turbed solution (a requirement which does not need to be 
considered when »>0, since orthogonality is assured by 
the presence of the factor P,,(cos @,.)). We have, therefore, 
(since «°(1) is normalized), 

‘ saiiae a 
a(R) = —(4nr/Z L | Pin" p*dp 

70 
- : 
t | ( out t B°u' up*dp (3) 
J R ' 

The final expressions for f,(1, 2) may be conveniently 

thrown into the form: 


For p<R, fxz=—(1/Z){A.(p)B,(R) 
+[C.(p)—D,(p) JE.(R H,(R)} 
For p> R, frx=—(1/Z){C,(p) LE. (R) 
—F,(R —G,(p R*-4 K R;} 
For n=0, we find 
Ao(x)=1 By(x) =0 


Do(x) =e** /2x— Ei(2x) 


F(x) 


Co(x) = 1/2x—x—log (2x) 
Eo(x) = —e-**(x+1)/x =(x—1)/x 
Go(x) =x 


Ho(x) =e"** |x +-2C—%8+(2C—1)/x 


+(1+1/x) log (2x)+(1—1/x)e**Ei(—2x)}. 
Ko(x) = Ho(x) — | (e~** Ei(2x) — 1)(14+1/x) 
+log (2x)(1—1/x)} 
C represents Euler’s constant. 
For n>0, we let n(n+1)/2=m; then 
A,,(x) =(n+1)x*+nx"*1, B,,(x) =1/2mx"*", 


C(x) =1/m+1/x+m /2 !x?+-m(m—1)/3 'x3 


+-m(m—1)(m—3) /4'x4 


+m(m—1)(m—3)(m—6)/5!x5+---. 


3 


D,,(x) = fe?*(— 1) "+14 1/x®— (m—1) /x 


+ (m—1)(m—3)/2'xt*—---}. 
E,,(x) = e-*7(—1)"*1(2n+1)(1+x) {1/x 
+m /x?+m(m—1)/2'v3+---}, 
F,,(x) = — }(2n4+1)(14+<x) {1/x—m/x* 
+m(m—1)/2'!x3—---+}, 


G,,. (x) =(n+(n4+1)x)B, (x) 


H,,(x)=0. K,,(x) =0. 
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Wave functions of several forms have been applied to the 
treatment to the ground state of Li. It is found that errors 
in the representation of the K shell have but little effect on 
the computed ionization energy, and that a satisfactory 
orbital for the Z electron can be constructed in a very 
simple form. The best wave function in which the relative 
positions of the K and L electrons are not taken into ac- 
count gives an ionization energy in error by 0.05 ev. To take 


INTRODUCTION 


HE lithium atom, a simple system with all 
the essential characteristics of the other 
more complicated and perhaps more interesting 
alkali metal atoms, has been the subject of several 
wave-mechanical investigations.' All these treat- 
ments have been characterized by the con- 
struction the function from 
one-electron orbitals, and may accordingly be 


of atomic wave 
considered as more or less satisfactory approxi- 
mations to the Fock procedure.* In each case the 
computed atomic and ionic energies are in error 
by some 2.0 ev, while the errors in ionization 
energy (obtained as the difference of the com- 
puted atomic and ionic energies) are of the order 
of 0.1 ev. Wilson has reduced this error to 
0.044 ev using for the K electrons the orbital 
e~*-69r and for the L electron 1.345re~°-°*5" —e-!-5”, 
That so simple a method of treatment should 
yield such satisfactory values of the ionization 
energy is a matter of considerable interest. The 
reliability of these values must suffer, however, 
from their having been obtained as the differ- 
ence of two quantities, each of which was much 
less accurately known; there is also no variation 
principle to assure us that the true value of the 
ionization energy is an upper limit monotonically 
approached by the calculated values as the wave 
functions improve. 
We have therefore 
amining the effect on the ionization energy of 


been interested in ex- 


* A preliminary notice of a portion of these results ap- 
peared in Phys. Rev. 47, 700 (1935). 

'E. B. Wilson, Jr., J. Chem. Phys. 1, 210 (1933) and 
references given there. 

2 V. Fock, Zeits. f. Physik 61, 126 (1930). 


account of inter-shell polarization effects the Hylleraas 
method has been extended, and a wave function has been 
constructed which depends on the relative positions of all 
the electrons. A total energy in error by 0.068 ev (0.034 
percent) has been computed, and the theoretical value of 
the ionization energy is fixed as 5.363 +0.007 ev against an 


observed 5.364 ev. 


various approximations in the atom and ion 
wave functions, in determining the theoretical 
ionization energy with considerable exactness, 
and in the of the 
energy and wave function of the complete three- 


improving determination 
electron system. Various steps in the passage 
from the Wilson function to the most flexible 
possible form readily suggest themselves. We 
have examined the following, which appear to 
be the most important: (a) improvement of the 
separate orbitals, for which Wilson used func- 
tions of somewhat limited flexibility, (b) im- 
provement of the of the K 
electrons by removal of the requirement that 
they be represented by orbital functions, and 


representation 


(c) complete abandonment of orbitals in the con- 
struction of the wave functions of the systems. 
In the last two stages of our work we have 
employed the method of Hylleraas, which has 
been eminently successful in the treatment of 
the normal and first excited states of He,*: * ° 
which have many features in common with the 
ground state of Li. In principle the extension to 
the more complicated system is quite simple, 
but the practical difficulties encountered are 
rather serious. A description of the way in which 
these were overcome is reserved for later sections 
of this paper; we proceed now to a presentation 
of the we have 


and discussion results which 


obtained. 


A. Hylleraas, Zeits. f. Physik 54, 347 (1929). 
A. Hylleraas and B. Undheim, Zeits. f. Physik 65, 
759 (1930). 

5A. S. Coolidge and H. M. James (see paper in this 
issue, p. 676). 


*E. J 
‘E. 
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IMPROVEMENT OF THE ORBITAL FUNCTIONS 


The wave functions hitherto used in treating 
Li have been of the form 


K(1)6(a@/o;) K’(1)6(8/o;) L(1)6(a/o;) 
Y=|K(2)é(a/o2) K’'(2)6(B/o2) L(2)6(a/o2)|, 
K(3)i(a/o3) K’(3)6(B/o3) L(3)6(a/a3) 


where K, K’ and L are the space parts of K and L 
orbitals, while 6(a@/o,), etc., are the spin functions 
introduced by Slater.® Actually AK(1) has always 
been taken to be identical with K’(1), as would 
be the case in the ion, but is not required in the 
atom. Expanding the determinant above, we find 
¥ =K’'(1)|K(2)L(3) —K(3)L(2)}6(8/01)8(a/o2) X 
5(a/o3)+two similar terms which differ by 
permutations which assign different spins to the 
electrons. Because of the different spin factors, 
no combination of two different terms can make 
any contribution to the matrix element of energy 
or unity corresponding to this function (so long 
as the Hamiltonian is assumed independent of 
the spin), so that these elements will be simply 
three times the contribution arising from any 
single term in combination with itself. Now this 
property can easily be seen to be a general one, 
possessed by any function having the symmetry 
appropriate to Li. We may write such a function 


Y= (1, 2, 3)6(B/01)5(a/o2)6(a/o3) 
+ (2, 3, 1)6(B/a2)6(a/a3)6(a/o;) 
+ 0(3, 1, 2)6(B/03)6(a/o1)6(a/ee), 


the only restriction upon g(1, 2, 3) being that it 
must be antisymmetrical in the coordinates of 
electrons 2 and 3. With this understanding, the 
whole discussion can be carried out on the basis 
of the function ¢g(1, 2, 3). 

For separate-orbital functions, we have 


g(1, 2, 3) =K’(1){ K(2)L(3) —K(3)L(2)}. 


> =~) * 


The Fock function is the best function of this 
form, and it is clear that Wilson’s function is not 
a very close approximation to it. His K orbital 
in particular was rather inflexible, and a treat- 
ment of He using such orbitals gives decidedly 
poorer results than the Hartree method, and, 
a fortiori, the Fock method. In the next section 
we shall deal with even more radical changes of 


* J. C. Slater, Phys. Rev. 34, 1293 (1929). 


the core function than improvement of the 
orbitals involved, and it will become apparent 
that no appreciable change in the ionization 
energy is to be expected from such a modification 
of Wilson’s function. We have therefore been 
interested only in the effect of increasing the 
flexibility of the L orbital, changes in which 
will be more strongly reflected in ionization 
energy. We find that with an L orbital of the form 
L(1) = cyrye 989571 + coe! 2571 + get 571 + ge! 7511 
improvement of just 0.002 ev over Wilson’s value 
could be obtained. An attempt to make a similar 
improvement in his function by replacing the 
term 7 ,e~°-65" by Cre" "1+ Core” 92"! +-Carje° 8"! 
gave no gain at all. Evidently Wilson’s Z orbital, 
though of very simple form, is about as good as 
any which might be used in functions of this 
type; of the residual error in his ionization 
energy only 0.002 ev can be attributed to lack 
of flexibility in this orbital. 


IMPROVEMENT IN INNER SHELL REPRESENTATION 


In the description of Li, the use of orbitals in 
connection with the KA electrons is clearly in- 
adequate, as it is in the treatment of the 
analogous He. The representation may be im- 
proved by similar methods, thus giving us a 
function of the separate-shell type, for which the 
general form is 


g(1, 2, 3)=KK(1, 2)L(3)—KK(1, 3)L(2), 


where KK(1, 2) is an unrestricted function of the 
relative coordinates of the two electrons, having 
the general nature of a wave function for Li*, 
but not necessarily symmetrical in 1 and 2. 

We have computed the energy of a number of 
functions of this form, both with the simple 
outer function L=re~®™’, and with the better 
one L=1.4re~°*r—e-!*", which is nearly the 
same as Wilson's. Table I contains the results, 
as well as a description of the core functions 
investigated. With the exception of E, the core 
functions will be recognized as of the Hylleraas 
type, modified as discussed in our paper on 
excited He. 

Function E is of a new type, which can be regarded as 
derived from a Hylleraas function by expanding the powers 
of rie as functions of r;, 72, and cos @;2, and then selecting 


a finite number of the resulting terms, with coefficients to 
be varied independently. In these terms there are dis- 
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Variation of ionization energy with inner shell 
functions. 


TABLE I. 


Energies in electron volts. 


INNER FUNCTION ALONI COMBINED WITH P COMBINED WITH Q 





Error in Error in 

Sym- Ion ioniza- ioniza- 
bol energy Error Atom tion Atom tion 
{ —192.89 +4.22 — 198.19 —0.060 — 198.27 +0.019 
B —195.49 +1.62 — 200.78 —0.077 — 200.83 —0.028 
( —195.53 +1.57 — 200.82 —0.080 — 200.85 —0.048 
D —196.30 +0.81 — 201.56 —0.095 
E —197.02 +0.09 —202.28 —0.094 —202.32 —0.057 
F — 197.05 +0.05 — 202.32 —0.093 

Ideal —197.10 0.00 (—202.37) (—0.093) (—202.39) (—0.056) 
The functions used are as follows (not normalized): 

Ae rites), 

Bae rite) (ey ea(ri tre) J. 

e e 2.6875 (ri +re 

D e-8ritrD) Fey +ee(ri tre) +earire tere tre) j. 

Ee 8¢si2t 912) ey +2812 +C35127 212 +4812? +5512 


+6512 COS O12+¢75127! COS? O12 ]. 
tr) eytea(ritre) +esriretca(r? +r2) 
+O5li2 +r iz? |. 
P rse 0 S63, 
Q 1. Arze~° 73 — e151, 
The coefficients c, were determined so as to minimize the ion energy. 
In constructing the atomic function, the same values were used. 
The symbols gi2 and siz in function E represent, respectively, the 


greater and the smaller of r: and re. 
The last row is an extrapolation, in which the error in the ion energy 


has been made to vanish. 


continuities of derivative when r;=rz which make finite 
contributions to the kinetic energy, but do not destroy 
their usefulness. We have applied such a function to the 
treatment of the ground state of He, and with the inclusion 
of nine terms have reduced the error in the computed 
energy to less than 0.02 ev. Since this error is some three 
times that given by Hylleraas’ six-term function and the 
computations are no easier these functions are not recom- 
mended as a substitute under such conditions. It was 
found, however, that when it was desired to add a third 
electron to the system and to neglect polarization terms 
with this electron (as in the present case) the treatment was 
more conveniently carried out with this than with the 
Hylleraas function; hence the present computation. This 
ceases to be true, however, when inter-shell polarization 
effects must be included. Since functions of this type seem 
to be restricted in usefulness and the methods of handling 
them in the three-electron problem are quite complicated 
we omit any further description of them. 


It is apparent from Table I that, for the given 
L orbitals, the ionization energy comes out 
almost independent of the representation of the 
inner shell, so long as the latter is not too bad, 
and that the small changes which do occur are 


similar with both orbitals. This would un- 


doubtedly hold for any good L orbital. An 
obvious corollary is that when L has been given 
the best form for use in conjunction with some 


AND 
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as: = 


reasonably good core function, this will be very 
nearly the best form for any other good repre- 
sentation of the inner shell. As an example, we 
note that the linear coefficient in function Q is 
essentially the same whether it minimizes the 
atomic energy with core C or with core E. Now, 
computation with core function C has shown 
that Wilson’s L function gives a result better 
than function Q by 0.004 ev, and that a further 
0.002 ev can be gained by still further im- 
proving the outer orbital. In view of the above 
discussion we may confidently expect that the 
values in the last column of Table I (at least 
the last few entries, in which we are primarily 
interested) would be likewise raised by 0.006 ev 
if Q were replaced by the best possible Z function. 

The core functions from A to F are pro- 
gressively better approximations to the true ion 
function. It is of interest to attempt an extrapola- 
tion to determine the energy which would be 
obtained by combining this function with various 
outer orbitals. The last row of Table I gives such 
estimates. In particular, if the best possible L 
function were used, the residual error would be 
very close to —0.050 ev, which is evidently the 
error inherent in the forced separation of the 
shells. 

In the examples so far given the core function 
and the ion function have been identical, and 
there seems to be no doubt as to the propriety 
of computing the ionization energy as the differ- 
ence of the computed energies of atom and ion. 
The best function for the atom core having a 
given form will, however, presumably differ from 
the best ion function of the same form; in fact, 
as has already been pointed out, the ion function 
has different symmetry properties from the best 
core function. Thus there arises the question as 
to how the ionization energy shall be determined. 
In order to investigate this point, we made 
computations differing from the F—P combina- 
tion of Table I in that, in the first case, the 
coefficients c, were allowed to take on values 
minimizing the atomic energy, and, in the second 
case, in that this adjustment was made after 
separate coefficients were assigned to the terms 
in 7; and fo, and to those in r;? and r,. In the first 
case the energy was reduced by 0.0003 ev, while 
the second charge permitted a further improve- 
ment of 0.0002 ev. We conclude that there is no 
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significant difference between the best ion func- 
tion and the best core function obtainable with a 
given series of the sort which we use. This is 
fortunate, for it relieves us of the necessity of 
determining the coefficients of the ionic function 
and introducing them explicitly as restrictions 
upon the variation carried out to give the corre- 
sponding atomic function. We need only solve 
the secular equations for atomic and _ ionic 
problems, and take the difference as the ioniza- 
tion energy. It also enables us to save labor by 
assuming that terms which in the ionic function 
must have the same coefficient may be treated as 
preserving this relation in the atom, thereby 
forming, for of variation, a single 
term. 


purposes 


PHYSICAL INTERPRETATION 


Qualitatively, these results are easy to under- 
stand. If we had complete flexibility in both KK 
and L, we should arrive, on minimizing the total 
energy, at a definite function of a rather interest- 
ing type, in which the relation between the L 
electron and the K shell taken as a whole is 
formally that between the Z and K electrons in 
excited He, particularly in the 1s2s*S state, 
according to the Fock approximation. In line 
with our discussion of the Fock function for the 
corresponding singlet state (no treatment of the 
triplet state being available), we may expect 
the defects of this function to be, speaking 
roughly, that it neglects the polarization of the 
core by the valence electron and the consequent 
contraction of the outer orbit, thereby yielding 
too high a total energy. We would also expect 
that the K shell function will be but slightly 
different from that of the free ion. Both the dis- 
tortion of the ion function and the polarization 
effects should be less pronounced in Li than in 
132s *S He, because the core electrons are more 
firmly bound. Now, in 1s2s'S He, we have 
found that the defects of the Fock function are 
responsible for an energy error of 0.078 ev; this 
should be less for the triplet state because the 
electrons are to some extent kept apart by the 
antisymmetry. The error in the total energy of 
the best separated-shell function (indicated by 
SS) for Li should then be a few hundredths of a 
volt only, and this will also be the error in the 
ionization energy if we disregard the difference 


STATE 
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between the correct ion function and the KK 
factor in our atomic function. 

As computed for a separated-shell function the 
energy is roughly separable into energy of the 
core, energy of the valence electron in the average 
field of the core, and exchange energy, of which 
the last two parts give the ionization energy. 
Starting with the function SS, we may replace 
the core factor by another which may contain 
large errors with respect to the mutual relations 
of the core electrons, but which has substantially 
the same average field. While the total energy 
will be seriously affected, the ionization energy 
will be influenced only through the variation of 
the exchange term, a second-order effect which 
might be of either sign. On the other hand, if 
the L orbital is varied the change in the ionization 
energy will arise principally from modification 
of the energy of the valence electron, with the 
change in exchange energy making a second-order 
contribution. Thus the energy is 
essentially a characteristic of the L orbital used, 


ionization 


with a minor dependence on the core function 
through the exchange energy. When functions 
permitting separation of the shells give ionization 
energies in error by less than 0.05 ev, as in 
Wilson’s work, this cannot be taken as evidence 
that the function is a good one; thus, it is clear 
from Table I that some combination of A and B 
could be found to give exactly the correct 
ionization energy. This effect is due to the com- 
pensating errors in the exchange energy which 
arise from a poor core function. To obtain a 
more accurate value for the ionization energy 
which is also reliable one must deal with the 
mutual polarization of the K and L shells. 
It is to this problem that the next section of this 
paper is devoted. 


INCLUSION OF POLARIZATION EFFECTS BETWEEN 
THE SHELLS 


For the computations of this section, we have 
extended our modification of the Hylleraas 
method, constructing new terms in our trial 
series by multiplication of the terms already 
found most useful by powers of 73 and fe; 
and/or by suitable negative powers of r; (assign- 
ing the number 3 to the L electron). The new 
terms (as well as all of those previously discussed 
except those of the function EZ) can be specified 
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as follows: We choose as variables p;= «ri, 
p23 = kfos, in which a possible additional degree of 
flexibility is gained by introducing the scale- 


constant x. Let 


CLR, m, n/y, 6, 6] 


L 4 


= (x? /41)?e "Pz 9P1) 93" po” pi", 
tlk, m, n/y, 6, 6 | 
= (x®/47)2e beste) nok 93” py" 


Our terms in the series for ¢(1, 2,3) then have 


the forms 


SLR, m, n, PJ = prs? SLR, m, n/y, 6, 6] 
glk, m,n, p]= pie’ SLR, m, ny, 6, 6] 
hk, m,n, P}= pis’ Rk, m, n/y, 6, 6) 

SC R*, m, n, p |= pos” CLR, m, n/y*, 6, 6] 


In the bulk of the work of this section, the con- 
stants had the values 6=3, y=0.65, y*=1.5, 
k=1. We postpone for the present a discussion 
of their selection, as well as of the choice of 
appropriate terms. The most general extension 
of the Hylleraas method would, of course, require 
provision for terms with powers of two or three 
interelectronic distances occurring simultane- 
ously, which could not be represented in the 
above notation. We have avoided this com- 
plication. We have also, in accordance with 
the results of the previous section, treated 
(glk, m,n, p|+e[k, n,m, p]} as a compound 
term with a single coefficient. [t will be seen that 
the terms g are just those occurring in the 
previously discussed functions of the separate- 
shell type. The new terms f and / allow for 
detailed treatment of the reactions between the 
outer electron and the inner one having the same 
and the opposite spin, respectively. 


It might appear that so long as one’s interest is restricted 
to the ionization energy it would be possible to ignore the 
terms g[k, m, n, p], >¥0, which permit an accurate treat- 
ment of the core. Thus one would regard the atomic func- 
tion as constructed from a core function, depending on r; 
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and re, and an L orbital depending on 73, 7:3, and 722, and 
would argue from Table I that the ionization energy should 
be practically unaffected by change in the core function. 
But it must be noted that we have no guarantee that results 
parallel to those in Table I would be obtained if the outer 
function were freely variable. For rj2, 7:3 and re3 cannot 
take on their magnitudes with complete independence, and 
from this arises the possibility that (in a very awkward 
manner, to be sure) some defects in the core function might 
be compensated for by adjustment of the inter-shell polar- 
ization terms. Thus with a core function poor enough to 
bring this effect into play the calculated ionization energy 
might appreciably exceed the actual value. It is accordingly 
necessary to carry out the work with core and atomic 
functions so good that such second-order terms can be 
safely neglected, and this seems to require a fairly accurate 
treatment of the effects on each other of the K electrons. 


As a basis for our work we have chosen the 
atomic function constructed from the Hylleraas 
core function F and the Wilson JL orbital, 
omitting, however, four of the twelve terms with 
an estimated error of 0.001 ev. To this function 
we have added nine new terms designed to take 
account of the polarization effects, obtaining a 
total energy in error by 0.068 ev. It was possible 
to make a fairly definite estimate of the limit 
attainable by the inclusion of indefinitely many 
such terms, including those with higher powers 
of ri3 and res; and combinations of such powers, 
according to the general method which we have 
previously explained.’ Such an estimate indicates 
a convergence limit in excess of the experimental 
value by 0.051+0.007 ev. Now, this is, within 
the limit of uncertainty, exactly the error which 
we ascribed to the defective core factor of the 
function before the polarization terms were in- 
cluded,—that is to say, it is the error of the best 
free ion function which can be constructed out 
of the terms which remain when we write down 
the terms in our best atomic function and then 
strike out all factors referring to 73, 713, OF fe3. 
It appears then that to a satisfactory degree of 
approximation, the core function retains its 
identity even in the correct, nonseparable atomic 
function; that our atomic function is so good an 
approximation that the errors due to imperfect 
representation of the inner shell and to incom- 
plete treatment of the polarization are additive. 
To the extent that this is true, we may speak of 
the ionization energy corresponding to our func- 


7H. M. James, A. S. Coolidge and R. D. Present, J. 
Chem. Phys. 4, 187 (1936). 














tion, even though it is of a nonseparable type. 


Our best calculated ionization energy is thus 
5.347 ev, and the extrapolated value 5.363 
+0.007 ev, as compared with the experimental 
5.364 ev. 

We have not determined the coefficients of our 
best 17-term function. We give here a com- 
paratively simple normalized function which is 


very nearly as good. 
g(1, 2, 3) =8.70800g[ 1000 ]+ 0.46443 { gf 1010 ] 
+¢[1100]} +3.74530{¢[1020]+¢f1200}} 


— 2.25446¢[1110]+4.66553¢[ 1001 } 





— 0.6037 2g[ 1002 ]—7.76625¢[0*000 | 
— 1.85902 { g[ 0*010 ]+ ¢[0*100 }} 
+ 1.57214f/[0001 ]+ 1.93299h[0001 }. 

It will be recalled that the complete wave func- 
tion is the sum of three such terms, differing by 
systematic permutations of the electron numbers, 
each multiplied by a different spin factor. 
Accordingly the above expression for (1, 2, 3) 
has been normalized to 1/3. 

The energy given by this function is as follows: 


Total energy, Free ion, Ionization 


| Ri jh, Ri ih, (ev) 
Experiment — 14.95784 — 14.56160 5.364 
Function — 14.95215 — 14.55776 5.339 
: Error 0.00569 0.00384 —0.025 
: Percent 0.038 0.026 0.47 


Of the total “‘polarization energy” of 0.05 ev 
about three-eighths comes from the pair of 
electrons with the same spin, and five-eighths 
from the other pair. This inequality is in agree- 
ment with the general rule that polarization is 
more important in singlet than in triplet states 
of two-electron systems. 


METHOD OF CALCULATION 


The calculation of the matrix components of energy and 
unity hinges upon the evaluation of the integrals 


I(abcdef /aBy) =| d Ve art Bret) 11279" 73° 19347 3)°r 07 . 


Consider first the case that d, e, f are each either 0 or —1. 
If e (say) is 0, then the integral is readily evaluated by 
taking as variables of integration rj, 72, 73, COS O12, COS 823, 
¢i3 (the angle between the plane containing r; and re and 
that containing rz; and re) and three unspecified angles 


giving the orientation of the whole system in space. The 
last four variables do not enter the integrand. We set 
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cos 923) and 1/rjg=Z,,5)2"g12 


1 /ro3 = = ,,523"g23 "'P, 


P,,(cos 0,2), where s;;, gj; mean, respectively, the smaller 
and the greater of 7;, 7;. Upon integration over the cosines, 


the terms for »>0 disappear, and we have, for example 


I(abc —10—1/aBy) = 649° { Wass. o41. 41(e, 8B, y) 

+ W, 19. «42, ofa, y, 8) + +2, atl, .4+1(8, a, ¥) 

t Wos2. 41 11(8, vy, a) + Wee. ase. oly, a, B) 
4-1 7, B, a)} 
where 


Wy, o. ala, B, =f. dx { dy [ dzx’e "6 Puc" 


I(abc — 100 aBy) =642°A, ,o(a) 
| Vore, 4118, y¥) + Vere, nurly, B)}, 


where 


A ;(a)= | dxx’e~™ 
: /0 
and 


fa n B, 4 )= {ay {" dz) me PUSN¢ 


a)Vm. »(B, y) = W;. m. ala, B, y) 
+W, j, n(B, a, ¥) 4 Won. n B, y, a). 


We note that A, 


If d=e=f=—1, it does not seem possible to find a sym- 
metrical method of evaluating the integral. We have ex- 
panded 1 /r;3 as before in terms of P,,(cos 613), which, by the 
theorem of biaxial harmonics, can be replaced by P,,(cos @;2) 
P,,(cos 693) plus terms which disappear upon integration 
over ¢}3. We thus find 
I(abe —1—1—1/aBy) =642°2,,(2n+1) 

X | Wasonte, 641, e—2n(a, B, y) + Wasense, o4 9 
+ Ws ,att,c 2n(B, a, y¥) + Woiense, c41, a—2n(B, y, a) 


+ Wesons2. a+1, b—2n(¥, + Wi2n+2. b+1. a-2n(y, B, a)} 


1, b—2n\Q, Y, 8) 





a, Bp 


The final formula is thus symmetrical, and the summation 
converges rapidly. 

Because of our not considering individual terms in 
which more than one interelectronic coordinate explicitly 
appear, we did not encounter the case that e>0. (This case 


could, however, be readily handled by extending the 


following method.) By means of the relation 
“ee 
93° = 12° +13" — 2rors cos 603 


it is possible to reduce any case of positive d by steps of two 
powers of riz until one reaches either d=0 or d= —1. The 
resulting expressions will have the same form as the integral 
I, except that some of them will now contain powers of 
COS 693 in the integrand. These can readily be treated by the 
methods indicated above. The same reduction applies to 
the case that f>0. It was found advantageous to work 
out the explicit formulas, in terms of the W's, giving 
I(abcdef/aBy) for every set of numbers d, e, f which we re- 
quired, as this resulted in considerable cancellation, and 
also improved the rapidity of convergence when all three 
numbers were odd and a summation had to be used. 

The first step in the calculation is the construction of 
V, and W, according to the 
By the 
systematic use of the latter, the work can be made self- 


tables of the integrals A, 
following definitions and recurrence relations. 


checking by verifying that two independent chains of 
calculations, in which the relations are applied in different 
order, lead to the same final value. 
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A; a) = I, e “x/dx, j “0. 
Vinn(a, B) = | fe aryMePUy"dxdy, m=0, m+n=-—1. 
r<y 
Wala, B, y= | | Je arySe Puy /e ¥*ohdxdyds, 
T<y<e 
f=0, f+g=—-1, f+eth=-—z2. 
A o(a) == | a; 
A ;(a) =(1 a)’A, (a), j=l 
’ 1 428 
Vo, -1(@, 8) = log — : 
a B 
° Fé n(a, 8)=(1 a)[A,(8) —A,(at 8) ], n—0; 
Vin, n(a, 8) =(1/a)[mVm_1, nla, 8B)—Amin(at8)], 
m=1, m+n=0; 
= (1 B)[Am in(a@ +8) +n Va. n—1(a@, 8)], 
m=0, m+n=0; 
Ricks hai " . 
= [mVm—1, n(a, 8)+nVim, n—1(a, 8) ], 
ats m=1, m+n=0. 


*Wo, g, ala, B, y) =(1/a)[Vo, a(8, ¥)— Vo, nlatB, v)], 


g=0, gth=-—l. 
Wy, g, hla, B, y=(1 a) LfWy 1, g, Ala@, B, y) 
— Veig, alatB, y)], f=1, f+eg=0, f+eth=—l; 


=(1/8)[gWy, 9-1, ala, B, v) 
+Vy+9, nlatB, y)— Vy, grnla, B+y)], 
f=0, f+g=0, f+gth=—1; 
=(1/y)[AW,, 9, rn-1(a, B, v) 
+Vy;, ginl(a, B+ry)], f=, 
1 
atb+y 


f+g=—-1, f+gt+th=—1; 


CfWy-1, 9, ala, B, vy) 


+gWy, 9-1, nla, B, y) +hW,y, 9, nila, B, v)], 
f=1, fte=0, f+eth=-}. 


By reversing these recurrence relations, the functions can be 
computed for negative m, g, and hk, within the limits for 
which they remain finite. It is well to avoid the relations 
marked,* as the computation of tables of all the W’s re- 
quired to construct all the 7's for a given choice of a, 8, y, 
can then be based upon preliminary tables of A(a+8+vy) 
and of V(a, B+y7), V(8, at+y), V(y, a+8), Viat+s, vy), 
V(B+y, a) and V(a+y, 8). 

Each term in our trial series consists of two parts, but 
because of the symmetry of the Hamiltonian operator, the 
matrix elements of energy and unity between any two such 
terms contain but two components, each taken twice. 
These components are computed separately. The elements 
of unity need no further discussion. Those of energy are 
found to fall into two classes, according as the interelec- 
tronic distances specifically involved in the combining 
functions refer to the same or to different pairs of electrons. 
The first class gives rise to the integral 


> ee ee , , , 7 
(29/6423) fd V { e~(" Pat BY pot’ pe) 5k’ 5m” 5M” 5?’ 
—(¢ +8 pb+ ; 27 4 
x Hie (a pat B pbt+ ¥ px oa" pe™ a, ™ oan? } =L +x«M;: 





AND 


COOLIDGE 


a> 
L=—(a’+6"+7")1(000, 000) +a(2k+2+ p)1(— 100, 000) 
—k(k+1+ p)I(—200, 000) +8(2m+2+ )) 
x 1(0—10, 000) — m(m+1+ p)I(0—20, 000) 
+-+(2n+2)I(00—1, 000) — n(n +1)J(00—2, 000) 
+p{|—(m+k+2+2p)1(000, 00 —2) —al(— 120, 00—2) 
+aJ(100,00 —2)+kJ(— 220,00 —2) —g7(2—10,00—2) 
+81(010, 00—2)+mJ(2—20, 00—2)} 


M = 1[000, —100]+J[000, 0—10]+J[000, 00—1] 
—3I[—100, 000]—37[0—10, 000]—37[00—1, 000). 


In these formulas, J[uvw, xyz] is an abbreviation for 


I(k+k’+u, m+m’'+2, n+n’'+w, x, y, 
p+p'+2/a+a’, B+8’, ¥+7']. 


The second class of combinations produces the integral 


’ a (a’ pat+B’ pht+v’ pc) 1 , , »” 
(2x9 /649°) fd V {e7 (@' Pat B’ wht’ pe) 5 RY jm? ym! yp) 
+ + ) a U 
« H {ee Pat Beoty pc pat pe pe" phe? 1 


The value of this integral is represented by the same formula 
as before, provided that J[uvw, xyz] is now understood to 


be an abbreviation for 


ILk+k’+u, m+m'+9, n+n'+w, p’+x, y, 
prez at+a’,B +p’, 77 7’ ]. 


These formulas are not symmetrical; the coefficients of 
the J's are determined solely by one of the combining 
functions. This is in practice rather an advantage, since it 
provides two different methods of calculating the com- 
ponents, and thereby checking the rather treacherous step 
of picking out the coefficients. A symmetrical formula 
could, if desired, be worked out by a method which we have 
discussed in a previous paper.’ The same paper contains a 
description of the process of setting up and solving the 
secular equation, and thereby determining the best func- 
tion which can be constructed from a given selection of 
terms. 


CHOICE OF PARAMETERS AND TERMS 


It would be desirable to have a form of trial function 
permitting independent variation of the screening con- 
stants for the inner and outer electrons, but this degree of 
flexibility could not be combined with that secured by the 
use of many terms with adjustable coefficients without 
tremendously multiplying the labor of computation. With 
our definition of the trial function, the screening param- 
eters are s;=«6é for the inner and so=xy, so*=«xy* for the 
outer orbital. The only practicable variation, that of x, 
changes all three in the same ratio. It is therefore important 
to start with suitable ratios of 5, y, and y*, and this can 
most easily be accomplished by making them equal to the 
values which 5), se, and sg* are expected to have, relying 
upon variation in « only for slight additional improvement. 
Accordingly, y and y* were fixed by rounding off Wilson's 
values, taking y =0.65, y*=1.5. To fix 6, we investigated 
the energy of the ion as given by a six-term function 
similar to F of Table I but with variable screening param- 
eter. Table II contains the results, which indicate that 


8H. M. James and A. S. Coolidge, J. Chem. Phys. 1, 
825 (1933). 











GROUND STA 
ase II. Computed energy of Li, as function of screeniy 
Ion EN \ I Ion ENERGY I ) 
Li Rij 
14.55844 0.043 3.0 14.55776 0.052 
85 14.55850 0.042 3.24 14.55452 0.096 


79 is the best value, and would give an energy error of 
0.041 ev, or 0.020 percent, to be compared with the 0.014 


error in Hylleraas’ treatment of He, in which the 





perce 


same six terms were used. However, when many terms 


the series are present, readjustment among the coefficients 


takes care of considerable variations in the screenit 


T 
4 





constant, which therefore need not be verv carefull 


chosen. We selected 6=3 for ease in computing. A smaller 
value would have improved our total energy by 0.011 ev 
but, in view of the results already presented, would not 
have affected our ionization energy appreciably. 

With 6, y, so fixed, the effect of taking « different 


from 1 can be anticipated 


and y* 
[he outer orbital will be directl 


nfluenced by the change in so and s9*, but since we have 


already chosen practically the best values for these 


constants, this will cause but a second-order change in the 
ionization energy. A poor choice of * could have 
been to a large degree corrected by varving «x. The effect on 


the core function of such a variation will be small, owin 


and 4 


g to 


readjustment of coefficients, and the residual effect will not 
be large enough sensibly to alter the ionization energy. We 
conclude that, so far as ionization energy is concerned, the 
chief usefulness of the variable « is to correct errors in 
choice of y and y*, and that in our case little gain could be 
expected. We have verified these relations by calculations 
with a single-term Z orbital, for which our y is not the best 
value of Se 
In selecting terms, we first set up a function for the ion 
using the terms occurring in Hylleraas’ function for the 
inalogous ground state of He. Upon multiplying these by 
the main term of the L orbital, re 
g(1000}, gf{1010}+ [1100], 
ef 1110], gf 1001}, and gf 1002) 
which gives an ionization energy error of 0.093 ev 
combination, Table I) 


we obtain the six 
g(1020 }+¢[ 1200" 
the best combination of 
F—Pp 
To these terms there correspond six 
g[0*010 ]+¢[0*100), etc 


from the minor term e “2 ” in 


terms 


“star’’ terms, g[0*000 ] , arising 
the Z orbital. From analogy 
with results with core function £, we should expect that 
inclusion of these terms would reduce the error to —0.056 


ev. [t was actually found that the first two star terms alone 
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sufhced to carry us within 0.005 ev of this value, and we di 
not go to the labor of yuting the matrix elements re 
qu red for the remainir g terms. We did find that inclus 
of the additional terms g{[ 2000] and g[1*000) produced a 
improvement of 0.006 ev, leaving an error of 0.055 ev 
Chis improvement is clearly due to perfecting the J 
orbital; it is exactly the same as was found with the crude 
core function C upon going trom orbital Q to the bes 
combination of six terms, which is a rather striking example 
of the degree to which corrections for different sorts o 
defects in the function are additive, and supports our cot 
clusion that the ionization energy error for the bes 
separate-shell function is —0.050 e\ 

lo take account of inter-shell polarization, we include 

the following additional terms 0001 1001. #7O*001 
h{ 0001), ATO011)+AT0101), kf 1001 ind Al O*001 Che 
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According to Herzfeld, the “chemical constant’’ of the electron gas is different in the high 


temperature region, where all actual measurem 
temperatures. Hence the value of the constant 
high temperatures from the well-known low 


present paper gives a calculation of the high 


quantities and a numerical estimation of these. 


1 


T is well known that the constant A in 
Richardson's formula! for the thermionic 
current 2 


t= AT*e¥'R! (1) 


W is the measured work function) should have 
under ideal conditions the value 120 amp. 
cm? °K. Under “‘ideal conditions’ we mean that: 

a) The reflection coefficient 7 for electrons 
should be zero, i.e., all electrons striking the 
surface should be absorbed by the same. 

b) The work function W should be constant. 
The work function is defined as the energy 
necessary to remove one electron from the top 
of the Fermi distribution and bring it to a point 
near the surface, but still to a distance from it 
which is large as compared with the lattice 
constant. The “‘top of the Fermi distribution”’ 
again means the energy of that electron which 
has the highest energy if all electrons are in the 
lowest possible state. This is the work function 
which is measured experimentally by the photo- 
electric effect if Fowler’s method? is used for 
the evaluation. 

(c) That the electron gas in free space should 
behave like an ideal gas. 

d) Finally, the surface must be uniform and 
thus free from patches and contaminations. 

This last condition and its effects have been 
carefully analyzed by J]. A. Becker.’ It seems 
that it can account in the case of the most 
careful experiments, and for low accelerating 
potentials for small deviations of A from 120 


‘For the history of Richardson’s discovery cf., e.g., W. 
Schottky, Handbuch der Experimentalphysik, Vol. XIII 
Leipzig, 1928), or A. L. Reimann, Thermionic Emission 
London, 1934). 

2 R. H. Fowler, Phys. Rev. 38, 45 (1931). 

3]. A. Becker, Rev. Mod. Phys. 7, 95 (1935). 


ents are carried out, from its value at very low 
A in Richardson's equation must be different at 
temperature value of 120 amp./cm? °K*. The 
temperature value in terms of thermodynamic 


only, and generally only for a lowering of this 
value. Under the conditions (a)—(d), (1) is a 
consequence of the principle of detailed balance 
and general quantum statistical laws.* 
Condition (c) is fulfilled for not too high 
temperatures® and will be assumed further on. 
If the reflection coefficient is different from zero, 
we should only expect an additional coefficient 
1—r in (1). Since 1—, is certainly in the order 
of magnitude of 1 and smaller than this, the re- 
flection can explain only a deviation of A from 
120 by such a factor.® In spite of this, the 
experimental values’ are in most cases rather far 
from 120. It must be emphasized, however, that 
these experimental values are still subject to 
comparatively large variations and in many 
cases hardly reproducible. They are most trust- 
worthy if the work function is measured photo- 
electrically also. It is only too common an 
experience that a new set of experiments changes 
the measured values of both A and W con- 
siderably, generally both in the same direction, 
which indicates that not so much the measure- 
ment of 7, as that of its temperature coefficient, 
is crucial. In spite of this, it is generally accepted’* 


*Cf. reference 1 and K. F. Herzfeld, Phys. Rev. 35, 
248 (1930), also J. H. Becker and W. H. Brattain, Phys. 
Rev. 45, 694 (1934). 

5Cf. M. v. Laue, Jahrbuch der Radioaktivitét und 
Elektronik, Vol. 15, 205, etc. (1918) and W. Schottky, 
reference 1. 

6 L. Nordheim, Zeits. f. Physik 46, 833 (1928); Proc. Roy. 
Soc. Al21, 626 (1928); R. H. Fowler and L. Nordheim, 
Proc. Roy. Soc. A119, 173 (1928) show that the reflection 
coefficient is very small in general, around 0.07. See also 
Langmuir and Jones, Phys. Rev. 31, 401 (1928). 

7Cf. A. L. Reimann, reference 1 and A. L. Hughes and 
L. A. DuBridge, Photoelectric Phenomena (New York, 1932). 
For Cb, H. B. Wahlin and L. O. Sordahl, Phys. Rev. 
45, 886 (1934), for Rh unpublished data of H. B. Wahlin. 

8’ | am much indebted to Professor L. A. DuBridge for a 
discussion on this point. Cf. also his monograph, Actualités 
Scientifiques et Industrielles, No. 268, Paris, 1935. 
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Pape I. Measured values of the constants in Richardsoi 


equation. 








Cs 162 1.81 | Ta 60 4.12 
Ba 60 2.11} Hf 14.5 3.53 
Ni 1,380 5.02 | Th 70 3.38 
Pd 60 4.98 | Re 200 5.1 
Pt 17,000 6.27 | Cu? 65 4.33 
Mo 55 4.15 Ag? 0.76 3.56 
W 60 4.53 | Au? 40 4.32 
Zr 330 4.12 | Cb 57 3.96 
Rh 35 4.81 


that the deviations of A from 120 are real and 
must be explained by the breakdown of assump- 
tion (b). The present paper will be devoted to a 
discussion of (b). 


2 


The consequences of the temperature variation 
of the work function have often been analyzed,° 
especially by Bridgman and Herzfeld. I shall 
follow the latter’s treatment quite closely. The 
reason this subject is taken up again, is that our 
somewhat improved knowledge of the metallic 
structure allows a somewhat closer analysis 
which shows that in addition to the effect, 
analyzed by Herzfeld (which is reproduced), 
there is another effect which tends to compensate 
for it in many cases. This is the reason, perhaps, 
that the measured A’s are not all higher than 120, 
as one may expect from Herzfeld’s paper. 

It has been emphasized by Becker and 
Brattain’® that the quantities measured in case 
of variable work function (and thus contained in 
Table I) are W* and A*, defined by 

W* = RT°*(d dT) |Ini/T®, 


In A*=Ini 7?+T7(d/dT) |In1/T? 2) 
=In 71/7T*+ W*/RT, 


where 7 is the thermionic current. 

The procedure to obtain theoretical expressions 
for W* and A* is to calculate first 7 and obtain 
then W* and A* by (2). 


3 


In order to calculate the current, one may 
first calculate the vapor pressure P of the 
electron gas at a point with zero electrostatic 


9P. W. Bridgman, Phys. Rev. 31, 90 (1928); K. F. 
Herzfeld, reference 4. 
10 Reference 4. 


potential by thermodynamics. If P is given, 
the number of electrons striking the surface can 
be calculated by elementary kinetic theory (m is 
the electronic mass) to be P/(2amkT)'. If r 
is the reflection coefficient, this gives a current 7 
per cm? of the metal. 


t= Pe(1—r)/(2amkT)'. 3 
The current from the metal is equal to this in 
equilibrium and is supposed to be equal also 
under the conditions of thermionic emission. 
Thus the problem reduces to the calculation of 
the vapor pressure P. 

Both the vapor pressure and the work func- 
tion!" depend, in addition to the temperature 7, 
on the (mechanical) pressure p, under which 
the metal is kept. We shall assume that the 
metal carries always only an infinitesimally small 
charge only. If we allow ¢ electrons per mole of 
metal to escape and keep the volume of the 
metal constant, the mechanical pressure will 
increase by ep,(v, 7). The energy necessary for 
the removal of ¢ electrons, when the metal is 
kept at constant volume, is «W,(v, 7). It seems 
to be simplest to apply the Clausius Clapeyron 
cycle by keeping the metal at constant volume. 
We then have” 


0 |ln P(v, T)/0T=W.(v, T 


which gives 


In P(v, T)=- - +(5/2) In 7+ 7 


j=(3/2) In (2rmk/h?)+1n 2. 


The well-known Sackur-Tetrode value of the 
chemical constant must be increased by In 2 
because of the spin. In order to determine ?p,, 
we can apply the usual cycle to a piece of metal 


1! For the definition of the work function see the pre- 
ceding section. The vapor pressure P is a function of the 
mechanical pressure p for ordinary vaporization also, but 
this effect is rarely considered. 

2 In this section, temperature 7 and volume 2 will be the 
independent variables. Thus temperature derivatives are 
always to be taken at constant volume, volume derivatives 
at constant temperature. 
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which is deprived of ¢€ electrons permole. The 


pressure is p+ep,, the energy E+e€W,., so that 
we have 
dE aW d 
4+ ¢ tpt+ep.=7 pb+ep.), 6 
rari od ol 
which gives in addition to the usual 
OE dv+p=Tdp al (7) 
valid for any substance, the required 
OW, dv+p.=Todp, a1 S) 


since (6) holds for every ¢«. From (8) follows 


I[GNER 


- > 


We shall now calculate W* and A* by (5), (3 


and (2). We have to remember, for this, that 
when the thermionic current is measured at 
different temperatures, not the volume, but the 
mechanical pressure of the metal is kept at 
constant value, namely zero. The temperature 
derivatives in (2) are meant to be taken 

constant pressure, therefore. Hence, if @ is th 
volume expansion coefficient, the molal vol- 


ume, we have 


ra] a 1 
W*= RT*( +a In 
a7 ar T? 


aw v, O 
.--— 
07 : . ‘ 
or W*=W,(v, T)+t0aTp 11 
Tr OW ,.(v, OW ,(v, 0) 4dr ; ; : 
+T| | . _ (9) as the comparison with (9) shows. For A’ 
. a On ae we have 
In (A* A(1—r)) = —W.(v, 0) ‘RT+W.(v, T)/RI 
.7 
+{ [W.(v, 7) -—W.(v, 0) dr /Rr2)+rap, R. (12 


Here A is the constant for the ideal case, 120 
amp. cm~? °K-*. The last term in (12) is that of 
Herzfeld, the other ones have been omitted by 
him, since the work function chiefly depends on 
the volume. Its contribution to A* is not neg- 
ligible, however. In general, p, is negative, the 


pressure decreases if an electron is removed. 
4 


For the evaluation of (12) we shall make 
assumptions very similar to those of Herzfeld. 
The total energy of the metal contains two parts. 
The first, (v) arises from the motion of the 
electrons, depending, therefore, on the volume 
only. The second arises from the motion of the 
nuclei. The first part is the most important for 
the actual value of the work function, its change 
by the removal of electrons from the metal gives 
practically all of this quantity. For the second 
part, an expression RTD(© 7) will be used. In 
this, © is a function of the volume and will be 
changed also by the removal of electrons by the 
amount e©, if « electrons per mole are removed. 


>We shall neglect the small specific heat of free electrons 
cf. A. L.. Reimann, reference 1 


The total energy of a mole of the metal, out of 
which ¢€ electrons have been removed and brought 
to a point near the surface but still sufficiently 
distant from it in order to make the image 
force negligible, is thus 
E(v, T) +eW,.(v, T) = Viv) +€W,, (2, 0 

+RTD(O T)+eRO.D'(0 T). (13 


These assumptions can be justified on the basis 
of the usual assumption of fast electronic and 
slow nuclear motion," but I shall not enlarge 
upon this subject. From (13) one readily obtains 

W.(v, T)=W,(v,0)+RO.D'(0O T), (14 


which gives with (12 


A* (-) (-) 
In =—)’ 
A(l-—r) T T 


=(0, T)D'(0,T)—(0,./0)D(0/T)+2ap,/R, 


T 
+0, D’ ‘: 
e 0 7 : 


In (A*/A(1—r))=—C,0,./RO+vap./R 
=—30, O+ap./R, (15 
i( Ric on. H. Pelzer and E. Wigner, Zeits. f. physik 


Chemie B15, 445 (1932); F. London, Zeits. f. Physik 74 
143 (1932 
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where C 3R is the spec ific heat of the metal 
at constant volume, it has assumed its classical 
value at the temperatures under consideration. 
The thermionic emission can be considered as 
reaction, and cal- 


a generalized dissociation 


culated, hence, according to the transition state 
The 


here, 


method for calculating reaction rates. 


result is identical with the one obtained 


but I preferred this method because it is the 
in these problems. As has 


more usual one 


been remarked by Herzfeld, (15) corresponds to 
the expression for the chemical constant of the 
the 


electron high while 


familiar 120 amp. °K? cm? is the low tempera- 


gas at temperatures, 
ture expression. This is quite similar to the 
at low temperatures the 
2rmk*8 /h*) while 


additional 


situation for Hoe, e.g.: 
chemical constant is (3 2) In 
at high there is an 
In 8x°kJ ‘h?. In contrast to He, the change of 


temperatures 
the ‘‘chemical constant” in our case is due to a 
change in the condensed state. Also the tempera- 
ture dependence does change in H» (from 7°”? to 
T?7?) while there is no such change in our case. 


5 
Next we shall try to estimate the quantities ©, 
and p, occurring in (15). It must be remarked at 
the outset that these two quantities are of a 
rather different nature: ©, is a function of the 
volume only, but p, depends, according to (9), 
on the temperature and can be reduced to the 


more basic quantities, dW,(v, 0) dv, ©, and 
d@,/dv by (14). One obtains thus 
In (A*/A(1—r))=—C,0,/RO 
— (va, ROW, (2, 0) dv 
—veal(d dv)(C,.0,. RO) (16) 


(9) is just the v derivative 
the 


since the integral in 


of the expression we have calculated in 


previous section. 


For the history of this method cf. E. Evans and M. 
Polanyi, Trans. Faraday Soc. 31, 875 (1935). The formula 


i= AT%>" RI 


used by J. A. Becker and W. H. Brattain (reference 4) has 
been derived in the papers to which they take reference only 
for the case when the work function is independent of 
temperature. Of course, it is always possible to represent 
the thermocurrent by a formula with variable w, but 
this w will not be equal to the work function. It follows 
furthermore from the third law that w can have no term 


- : 
linear in T at low temperatures (cf. reference 9). 
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EQUATION 


It would be rather difficult to estimate the last 


expression in (16) since it contains d@, dv. If 
one assumes however, that 7(d 07)(C.©, R®) is 
of the order of magnitude of C,.©, R© itself 

which seems reasonable, the last term can be 
neglected altogether, because al 10 Since 


at the absolute 


W,.(v, 0 


zero, to remove one electron, dW,(v, 0) dx 


is the energy necessary 
is the 
negative increase in pressure, due to the removal 
of one electron. Since the electrons do the binding 
in the metal, the removal of one will loosen the 
binding and the actual pressure decrease will not 
be quite as great as ¢ alculated on the free electron 
hypothesis. For alkalis, one can obtain an esti- 
mate of this the 
derived for the calculation of the work function 
Our W, (2, 
what is denoted there by ¢. Eq. (7) 


magnitude from formulas, 
0) is in the present approximation 
reference 16 


reads in our notation: 


—0.58/(r,+5.1)—0.197, (r,+5.1)?—eD). (17) 
The work function is expressed here in volts, r, 
in Bohr’s units (0.528-10-* cm) and is defined by 
4ar3 3=v9 L=atomic volume. D is the mo- 
mentum of the double layer on the surface and 


will be neglected hereafter. From (17) we obtain 


dv=(r,/3)0W,/dr.= 


v0 W, 
—4.5(rdEo/dr,—7.4/7r.2—1.18/r, 


+0.587, (r,+5.1)?+0.38r,2, (r,+5.1)*). 18 


In order to eliminate dE,» dr., we set the deriva- 
tive of (8) reference 16 equal to zero. This 
expresses that the total energy is a minimum 


for r.. 


—rdkEy dr,+4.42, r.2+0.28 4, 


0.58r,/ (r,+5.1)7=0. 19 
This gives 
vIW,(v, 0) / dv 
= 4.5(3/r.?+0.9/r, —0.387,?/(r,+5.1)8 20) 


which is for r,=4 about 1.5 volts. The second 
(16) a= 21.10 


— 3.6. 


term in thus becomes with 


about 


]. Bardeen and E. Wigner, Phys. Rev. 48, 84 (1935 
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The quantity ©, © is the percent change 
of the vibration frequencies for one percent 
change in the number of electrons. It is, there- 
fore, half the percent change of the restoring 
force F, if one percent of the electrons is removed. 

In order to obtain an estimate for this, we 
may consider first the whole restoring force if 
only one ion is displaced. This is the force which 
acts on the displaced ion and originates from the 
electron cloud which is itself distorted by the dis- 
placement of the ion. 

It is difficult, in general, to calculate the dis- 
tortion of the electron clouds. Naturally, the 
inner electrons can be considered to be rigidly 
attached to the ion. We may assume that the 
valence electrons have an even distribution 
throughout the lattice and show wriggles in the 
neighborhood of the ion only. Then, the total 
charge distribution of the valence electron can be 
considered to contain two parts: p;, the varia- 
tions of which are appreciable only in distances 
comparable with the lattice constant, and p»2 
which contains the small wriggles only. It seems 
reasonable to assume that the ion carries the 
wriggles with itself but as long as only one ion 
is displaced, the slowly varying part of the charge 
will not be much affected. If the ion was displaced 
by x in the X direction, it will be acted upon by 
a force 

Fx =xZedE,/ dx, (21) 
where £ is the field due to the slowly varying 
part p, of the electronic charge and Ze the ionic 
charge. We have, by Poisson’s equation div E 
=47p, and for a cubic lattice, this gives because 
of dE,/dx=0E,/dy=d0E,/dz 

F=4rp,Ze/3. (22) 
The percent change of this,!? for a one percent 
decrease of the number of electrons is simply 
— pir Pim Where py, is the slowly varying part of 
the density of the highest energy electrons, at 
the point where the ion is; p;,, is the mean value 
of the same quantity for all valence electrons. 
The first term of (16) is therefore, approximately 


(3 2)( pit Pim). (23) 


?? 


'7Qne can calculate, from (22), the characteristic 


temperature of the metal, and it comes out in the right 
order of magnitude. A more satisfactory calculation of the 
vibrational frequencies of a metal has been given lately 
by K. Fuchs, Proc. Roy. Soc. A153, 622 (1936). 
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This will be comparatively large (~1.5), if the 


last valence electron is an s electron, small 
otherwise. If the highest energy electron happens 
to be such an electron in the lattice, for which 
pi: is zero, the whole first term will be negative 


also, because of the approximations we made. 


6 


It is hardly necessary to mention that while 
the considerations of sections 2 and 3 are, (12) in- 
clusive, strictly based on thermodynamics and 
those of 4 also should be correct within a very 
small error, the considerations of section 5 are 
very crude and will not give, in general, more 
than the sign of the effects and their order of 
magnitude. What can be claimed safely, is only 

(a) that it is purely accidental if the constant 
A* in Richardson’s equation (defined by (2) 
has the value 120 amp.,/°K? cm? sec. 

(b) that the deviation in the In of A* should 
be of the order of magnitude 1-3, rather negative 
than positive. It should be negative always when 
the last valence electron is not essentially an 
electron in the lattice. 


It should be emphasized once more that rather 
large deviations in the constant A* can be caused 
by a patchy character of the surface, as con- 
sidered by Becker and Rojansky.* These patches 
can consist in surface contamination or also in 
the polycrystalline character of the surface, 
since different crystal planes will have different 
work functions. It would be desirable, from this 
point of view, to measure the work function and 
A* on definite crystal planes, i.e., to use single 
crystals for the experiments.'* 

In some cases, zero energy outside of the metal! 
may correspond to a forbidden region inside. 
This would have as consequence an unusually 
large reflection coefficient.'!* The effect of a 
n 


=. 


reflection coefficient has also been omitted 
this paper. 

It is a pleasure to express my gratitude to 
Professor Herzfeld for his valuable discussion and 


criticism of this paper. 


i$ Professor H. B. Wahlin has pointed out to me that 
part of the effect of the continued heating and flashing ot 
the samples is perhaps due to recrystallization. Cf. H. B 
Wahlin and J. A. Reynolds, Phys. Rev. 48, 751 (1935). 

19 P. M. Morse, Phys. Rev. 35, 1310 (1930). 
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1 


N his paper, A New Relativity, Page! introduces 

reference systems of particle observers where 
the comparison of times and the measurement 
of distances are based solely on geometrical 
optics. Hence, by means of optical experiments, 
two such systems will be indistinguishable if 
suitable coordinates (x, y, 2, f) can be introduced 
for the observers of each system such that light 
travels along straight lines with constant speed 
and the space coordinates (x, y, s) are Cartesian 
coordinates of a Euclidean space. These proper- 
ties are preserved by Lorentz transformations. 
Page considered the question of the existence of 
such reference systems which are not Lorentz 
transforms of one another. He answered the 
question in the affirmative by constructing sys- 
tems with elements having constant acceleration 
with respect to an original inertial system. He 
also proposed the problem of investigating other 
types of such systems and the present paper 
contains the determination of all these. It is 
further shown that the assumption of constant 
light velocity plus the existence of one reference 
light paths implies the 


system with straight 


straightness of light paths in all such systems. 


2 


Consider two reference systems S and S’ of 
the above type. A correspondence is obtained 
between S and S’ as follows. To each particle 
observer P and time ¢ in S we associate the 
incident particle observer P’ at the time ¢’ in S’. 
We consider the four dimensional space-time 
manifolds (x, y, 2, t) and (x’, y’, 2’, t’) describing 

and 5S’. 


spondence may then be written 


the observers in S The above corre- 


x’ =9,(x, y, 5, t), 2’ =@3(x, y, 5, 2), 
(T) 
, ; 
VY =O2(X, Vv, 5,¢), £ =4(X, v, 5, bt). 


We introduce the vector notation x=(x, y, z, ¢) 


=(X1, X2, X3, x3). The scalar product a-b of two 
vectors a and b is defined as 


Leigh Page, A New Relativity, Phys. Rev. 49, 254 (1936). 
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a-b=a,),+ 12027-4303 — C04) +5 


The light paths in S are of the form x=a+ Ab 
where a is a vector from the origin to a point on 
the path, b is an ’ vector, i.e., b-b=0, 
and X is a real parameter. Our problem is to 
under the hy pothesis 
of constant light velocity, ie., all light paths 
a+b, b-b=0 in S are carried into isotropic 
curves for which 


‘isotre pic , 


completely characterize (7 


curves in S’, i.e., 
dx-dx=dx*?+dy?+dz*—cdt??=0 


We may reduce the problem to known results 
by the theorem: 

If all light paths x=a+ Xb in S are carried 
into isotropic curves by (T), then (T) 
Ix?+dy*+dz?*—cCd?=0 in- 


leaves the 
differential equation 
variant. 

For since (7) carries a+b into an isotropic 
curve for all b-b=0, on substitution in (TJ) we 


see that 


3 ‘ 00, 4 00; 
) p b. | Fi > b, =) 1 
l Ox, 


1\j=10xX 


But dx is an isotropic vector if x(\) is an isotropi 
curve. Hence, replacing 0; in (1) by dx; we 
obtain dx’-dx’=0, or the transform x’(\ 
isotropic curve X(A) is isotropic. 

The problem is thus reduced to that of de- 


ot any 


transformations which 


If we allow 


termining all preserve 
dx*+dy*+dz*—cCdt =0. 
nates to vary over a complex domain, ¢ may be 
replaced by (—1)*(r « 
that of determining the transformations which 


the coordi- 
and the problem becomes 


leave dx*+dy?+dz*+dr*?=0 invariant. All such 
transformations were given by Lie? who showed 
that they are all conformal transformations and 
form a group depending on 15 parameters. He 
further showed that this group is generated by 
the orthogonal transformations which carry the 
expression (x —xX,)?+(x—Xe2)*+(x—4x3)?+(x—xy 

into a constant multiple of itself, and the trans- 


2S. Lie, Theorie der Transformationsgruppen (Leipzig, 
1893 
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formations by reciprocal radii. The latter have, 
of course, singular points which are carried to 
infinity. In the following paragraph we consider 
properties of the analog of these reciprocal radii 


transformations. 


3 
The analog of the transformation by reciprocal 
radii is 
x’ = P,(x)=x/(x-x). 2) 
This transformation is valid for all elements 


(x, v, s, 2) which are not on a light path through 
the direction 


the origin. The effect of (2) on 


dx dd of a tangent to a curve x=X(A) is given by 
dx’ ~( x ) 
dy dr\x-x 
1 r ax dx\ 7 
md (5 -E)—— 2Ee B 
(x-x)? dn dn/ . 


is a light path not 


. (3) 





In particular, if x=x(\ 
through the origin we have x(A) =a-+Ab, b-b=0. 
Hence dx d\=) and 


= [ (x-x}b—2x(x-b) | 


= [(a-a)b—2a(a-b) ]. 
(x-x)* 


Since the vector factor on the right is indepen- 
dent of A the the 
constant and we have the theorem: 

The image of a light path is a straight line. 

The speed of an element described by the 
four-dimensional curve x=x(A) is ¢ if and only 
if (dx dd)-(dx dd) =0. That the velocity of light 
is constant follows simply. However it is con- 
venient to calculate (dx’, dd) - (dx’ dd). We obtain 


direction of tangent is a 


dx’ dx’ 1 dx dx\ 7 
ee Sees IB canes | x» ~2x(x: ) 
dy dy (x-x)4 dx dn/ - 


4 


r dx dx 1 (dx dx 
[ex ~2x(x: )|- ( . ). (4) 
dx dx x-x)*\dd dX 


As an immediate consequence of (4) we conclude 





not only that the velocity of light is constant and 


that isotropic lines are preserved but the validity 


ROM 





ZORN 


AND M 


of the theorem: 

The quotient of two physical intervals (dx-dx 
(d\x-d x) 1s invariant. 
that the 
gether with the Lorenz transformations generat 


\ e shall de 


scribe this result more explicitly in the next 


Lie asserts transformations Py, 1 


all the required transformations. 


paragraph. 
+ 


In this paragraph we give the results of Lic 
in a form which is directly applicable. We com- 
bine Py with the translations 

aS 
x=— = P x), 
(x—a)-(x—a) 


which are valid for all elements not on a light 
path from a, and obtain 


C.(x) = PP (x) 


le) Ge) 


By a simple calculation it is seen that 





zZ-~Giz-s 
C,(x)= ; 
1+(a-a)(x-x)—2(a-x) 


Hence Cy is the identity. The singular points for 
these transformations, i.e., the points for which 
the denominator is zero, form a hyperplane if 
a-a=(), otherwise a quadratic surface. It can 
be shown that these singular manifolds are dif- 
ferent if a and b are distinct. Let us suppose that 
x’ = L(x) is a Lorentz transformation or propor- 
tional to one. It is known that: (1) L(x) has no 
singular points, (2) straight lines are preserved, 
(3) x-x is preserved except for a constant factor. 
Hence all transformations LC, preserve the differ- 
ential equation dx-dx=dx*+dy?+d2*—cCdf=0. 
Furthermore we have the following theorem: 

If L1C,=L2C, then L,;=L2 and C,=C, 

For the singular manifolds of ZC, and LeoC 
must be identical. Since Z; and Ly» are regula 
this implies that C,=C, and hence L,;=L»2. As 
a consequence we see that the transformations 
LC, form a fifteen parameter manifold, hence 
they contain all the required transformations 
sufficiently near to identity. 
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LETTERS TO 


Potential Gradients in Direct-Current Metal Arcs Deter- 
mined from Stark Effect Studies 


A recent survey of the theory of operation of direct- 


current arcs in air between low melting point metal 


electrodes indicates a rather wide-spread adherence to 


the belief that the electron supply at the cathode is 


furnished by some sort of field current emission mechanism 
In support of this view and adding to the confusion alread 

there is quoted a conclusion drawn 
Stark 


existing in the field 


by Nagaoka and Suguira in their paper on the 


effect broadening of metal lines in arcs.! This conclusion 


s th 


at from the observed broadening of lines, one can 
infer the existence of electrical field gradients of hundreds 


the arc 


of thousands of volts per cm in and especiall) 
near the cathode. In quoting these conclusions in recent 
vears, workers on the mechanism of arc discharge fail to 


ind Suguira was done in 
1924), so that the 


interpretations of their results require restatement 


note that the work of Nagaoka 


the early days of Stark effect theory 


} 


eari\ 


in the light of present day knowledge. Failure to note 


these changes leads to obviously contradictory and absurd 
conclusions. No doubt exists that the phenomena observed 
by Nagaoka and Suguira and interpreted by them as a 


Stark effect is actually a Stark effect. Svmmetrical broad 
ening is observed in just those lines where the first-order 
Stark effect would be expected while similar conclusions 
apply to the unsymmetrical broadening ascribable to the 
second-order Stark effect. 

The difficult, 


Nagaoka and Suguira does not lie 


conclusions of 


introduced by the early 
in their ascribing the 
broadening te a Stark effect. The difficulty lies only in thei 
considering the electrical fields producing the Stark effects 
to be due to /inear potential gradients along the arc column 
amounting in some cases near the cathode to hundreds of 
kv per cm. Such gradients while perhaps welcomed by 
adherents to a field emission theory of the arc mechanism 


are completely at variance with the other data. In many 


cases observed by Nagaoka the broadening of the lines 
extends nearly along the whole length of the arc, so that 
the first-order broadening 
The 


existence of such gradients in the arc is completely nega- 


well in the center of the arc in 


the fields calculated are in the tens of kv per cm. 


tived by probe measurements of Nagaoka and Suguira 


and many others in these and similar arcs. Again an 
integration of the measured gradients over the length of 
their arcs would lead to potentials across the arcs hundreds 
of times greater than the applied potentials. 

It is obvious that the explanation of this apparent dis 
crepancy lies in the mistake of ascribing the Stark effect 
to linear potential gradients down the arc. The absence of 
any structure in the broadening of the lines indicates at 
once the absence of any such clear cut field gradients. The 
fields responsible for the type of Stark broadening ob- 
served can today be entirely ascribed to the ionic force 
fields existing in the intensely ionized plasma of the arcs 
The fields of 180 kv/cm correspond to an ion density of 


about 2 percent. The variations in broadening down the 
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are length observed | Nagaoka and Suguira represet 


obvious and expected char n 


path It is possible that valuable data as to the ion densities 
in arcs could be drawn from such observations 
Dunningtor 
roadening of the Mg, Cd 


Zn lines in sparks observed 10 reakdowt 


explanation was invoked by Lawrence and 


in their explanation of the | 


second after | 


In this case the fields were 10° volts cm and the densit 


of the ionization was estimated as of the order of one-third 
of the molecules present 
LEONARD B. Logs 
W. A. HILLEBRAND 
H. E. Waitt 
R. N. VARNEY 
F.C. MILLEt 
I versit f < r 
Berkeley, Califor 
pril 1 1936 
Nagaoka and Suguira, Jap. J. Phys 3, 45 (1924 
Lawrence and Dunnington, Phys. Rev. 35, 405 (1930 


Note on 6 Quartz 


Harold 


form ot cr 


In a 


announced 


recent issue of this journal Osterberg 


the discover ol a new stalline 


quartz, called 6 quartz, which he found to be not piezo- 


electric at temperatures below 89.6°K. We desire to call 


attention to papers by L. Balamuth? and F. Rose’ in which 


these writers describe the composite piezoelectric oscillator 
method for measuring the elastic constants of solids at low 
and } 


excited, respectivel , to 


temperatures. It is there stated that X cut quartz 


cylinders can be piezoelectricall, 


longitudinal and torsional vibration at a temperature as 


low as 78°K 

L. BALAMUTH 
F. Rost 
S. L. Ourmpy 

Col Un 

\pr &. 1936 

Octerber P Rev. 49. 55 1936 

Bal th, PI Rev. 45, 715 (1934 

Rose, PI Rev. 49, 50 (1936 

An Acknowledgment 
Both L. H. Thomas of Ohio State University and H. P 


Robertson of Princeton have drawn our attention to a 


paper by E. Cunningham! in which the author derives the 


transformations for the components of E and H under an 


inversion as given in our recent paper.2 Cunningham, 


however, does not seem to have realized that the inversion 
relative acceleration, 


transiormation represents constant 


nor does he discuss the force equation. We regret that 

we were unaware of Cunningham's work at the time we 

sent our paper to the Physical Review. 
LEIGH 
N. I. 


PAG! 


Apams, JR 
Vale University, 
New Haven, Connecticut, 
April 13, 1936 


London Math. Soc. 8, 77 (1910 


Cunningham, Prox ' 
Phys. Rev. 49, 467 (1936 


2? Page and Adams, Jr., 
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Enrichment of Nitrogen in the Isotope N! 


In a previous letter! we described briefly the production 
and analysis of samples of methane which were enriched in 


the carbon isotope, C™, by diffusion in a Hertz apparatus 


Further diffusions with methane have given a considerably 


greater enrichment than that reported. Two samples 


showed an intensity ratio of 3.3 to 1 for the bands due to 
CeCe and CC 


atomic percent tor ¢ 


corresponding to an abundance of 13.1 
The head due to CC! 


and its intensity 


could also 
be measured on these plates,” was con- 
with that of C#®C 


by one of us to photograph the C! 


sistent his methane has been used 
N" bands under high 
dispersion, and from these a value for the mass ratio of 
the carbon isotopes has been obtained 

The method has now been applied to nitrogen gas, with 


satisfactory results. Abundance measurements were made 
by observations on the isotope heads accompanving the 
1,0 and 2,0 bands of the second positive system of Ne. 
According to Urey and Murphy,® N® is present in ordinary 
nitrogen to one part in 350, so that these heads due to 


NUN 


heads 


should have an intensity of 1/175 that of the main 
NUN" A sample of 300 cc of 
p=6 mm) from the first diffusion gave an intensity ratio 
of 1, 16.5 
of N®, 
10. With this gas it should be possible to measure the 
bands due to C2N®, 
the mass ratio of the nitrogen isotopes 


DEAN E. 


due to nitrogen 


This means an abundance of 3 atomic percent 


an enrichment by e factor of slightly greater than 
and so to obtain an accurate value ol 


WOOLDRIDG! 
Norman Bridge Laboratory of P 
California Instit { n 
Pasadena, Californ 


A. JENKINS 
Department « 
University « 


serke 


Wooldridge and 
2 Due to an error in the rey f 
etter, the triple dot « gnating tl ad was conside 
The other two | correctly marked, and the 
to C#C12, CLC, and C3C, should to the eye appear eq y 
Jenkins and Wooldridge, paper given at the Washingtor 
f the American Physic: i \pril 30, 1936 
G. Herzberg, Zeits. f hemie B9, 43 (1930 
H. C. Urey and G. M Phys. Rev. 41, 141 
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eads were 


19%) 


A Spectrographic Method of Measuring Carbon 
Dioxide Concentration 


Chemical methods of estimating the CO» conc entration 


in a mixture of gases, such as in air, require very accurate 


temperature control and are slow in action. For example, 


in a continuous flow method! used in this laboratory the 
temperature of the electrolytic cell must be held constant 
to better than 0.01°C and about 30 minutes are required 
to obtain an individual reading. The physical method of 
observing the heat conductivity of the gas also requires 
an extremely accurate temperature control. This latter 
method measures a physical property of all gases, thus 
when used it is necessary to know the other gases present 
and make the proper corrections for them. Recentl, 
November 21, 1935) a spectrographic method was tested 
in this laboratory and found to be very successful. 

This method makes use of the high opacity of CO» gas 


to radiation! n the region of one of its fundamental absorp- 


»>DITOR 


tion bands. The spectral region used is a band from 4.2 t 
4.34, so chosen because water vapor and practically 

of the other gases or vapors commonly found in air are 
this radiation. Radiation of this 


transparent to wave 


length is isolated from a suitable source (Nernst filament) 
by means of an infrared spectrograph and allowed to fall 
upon a vacuum thermocouple after passing through 

known optical path length of air contained in an absorption 
tube. The apparatus is calibrated by filling the absorption 
CO known 


concentration of COs and the galvanometer deflections 


tube with air free from then with air of a 


For the range of concentrations from zero to 
0.10 percent CO 


With this arrangement the sensitivity of the apparatus is 


observed 


an 86 cm optical path length is used 


deflection 
COs». 


such that a one-millimeter represents 
Shorter 


measurement ol! greater concentrations. 


one-ten 


thousandth of one percent path lengths 


enable the Iwo 
millimeters of pure CO, at N.T.P. transmit only 22 percent 
of the radiation 4.2 to 4.3,y. 

The particular advantages of this spectrographic method 
The sensitivity is as that of the best 


are: (1) as great 


chemical method, (2) it gives instantaneous measurements 
galvanometer response 5 seconds), (3) the measurements 
are independent oi the humidity of the air and also inde- 
pendent of the presence of practically all of the other gases 
or vapors commonly found in air, (4) for ordinary room 
temperatures and pressures the apparatus has very small 
temperature and pressure corrections (3°C or 7 mm pres- 
sure change cause about 1 percent change in sensitivity), 
and (5) it can be recalibrated in a few seconds, which 
enables the use of full sensitivity. The same arrangement 
can be used to observe small concentrations of other gases 
by selecting radiation of appropriate wave-lengths 
E. D. MCALISTER 

Dis m of Radiation 
Smithsonian Institution, 

Washington, D. ¢ 

April 15, 1936 
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Focusing Effect of Double Crystal Spectrometer 


The focusing effect of the double crystal spectrometer 
in the antiparallel position recently discussed by Parratt! 
was developed in detail by the author? and later used by 
A. H. Barnes.’ 

If A@ is the difference in Bragg angle between the x-ray 
which passes through the axis of the second crystal and 
the two rays will intersect at a 


any other x-ray, then 


distance from the axis equal to sin A@/sin 3A@ (approxi- 
mately one-third) the distance from the target to the axis. 
In addition to the focus “fixed in space” there is also a 


focus “‘moving with the second crystal’’ which is three 
times as far from the axis, i.e., at a distance equal to the 
x-ray path from the target to the axis. 

The advantage of placing slits or ionization chamber 
windows at these foci was pointed out by the author. 

Roy C. SPENCER 
The Brace Laboratory of Physics 
University of Nebraska, 
Lincoln, Nebraska, 
April 18, 1936 
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